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Abstract 


Electromagnetic  materials  characterization  is  important  in  the  design  of  systems 
that  interact  with  electromagnetic  waves.  Determining  the  constitutive  parameters  of  a 
material  is  a  vast  area  of  research  and  practice.  For  this  paper,  discussion  will  focus  on  a 
destructive  method  using  waveguides  in  the  frequency  range  of  6-18  GHz.  Traditional 
methods  to  perform  similar  measurements  include  using  coaxial  cable,  stripline,  focus 
beam  and  rectangular  waveguides.  This  work  will  use  Double  Ridged  Waveguide  (DRWG) 
to  compare  to  these  other  methods  and  will  discuss  the  attributes  and  drawbacks  of  this 
new  approach.  The  most  similar  method  utilizes  rectangular  waveguide,  so  the  primary 
focus  will  be  on  comparisons  this  method  with  DRWG.  The  significant  advantage  to 
using  DRWG  is  the  increase  in  available  measurement  bandwidth.  The  challenges  include 
sample  fabrication  and  increased  mathematical  difficulty  in  finding  the  cutoff  frequency  for 
DRWG.  These  challenges  are  addressed  and  measurement  results  are  examined. 
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CONSTITUTIVE  PARAMETER  MEASUREMENT  USING  DOUBLE  RIDGE 


WAVEGUIDE 


I.  Introduction 

Measurement  of  a  material’s  constitutive  parameters  are  important  to  the  design 
and  maintenance  of  any  system  that  will  interact  with  electromagnetic  waves. 
Determining  these  parameters  can  be  especially  useful  in  Low  Observable  (LO)  aircraft 
platforms  for  diagnosis  and  health  status  of  the  Radar  Absorbing  Material  (RAM) 
treatments.  The  constitutive  parameters  can  vary  with  frequency,  thus  providing  a 
challenge  for  measurements.  The  measurement  techniques  used  need  to  be  selected  in 
accordance  with  the  material  and  scenario  in  which  the  measurement  is  taking  place.  The 
technique  explored  in  this  thesis  will  utilize  a  Network  Analyzer  (NWA)  and  waveguides. 
Traditionally,  rectangular  waveguide  is  used  since  its  electromagnetic  properties  are 
well  known  and  the  mathematical  model  is  well  understood;  however,  the  measurement 
bandwidth  is  narrow.  To  gain  a  larger  bandwidth,  this  thesis  will  explore  the  use  of  Double 
Ridged  Waveguide  (DRWG)  in  place  of  rectangular  waveguide.  The  scope  and  limitations 
of  using  DRWG  will  be  compared  with  other  methods  of  measurement. 

1.1  Problem  Statement 

While  the  NWA  is  capable  of  recording  measurements  across  a  large  frequency  range, 
the  transmission  line  can  limit  the  available  bandwidth.  For  the  well  known  measurement 
method  of  using  rectangular  waveguides,  the  waveguide  itself  is  the  limiting  component. 
In  order  to  acquire  a  broad  band  of  measurements  with  rectangular  waveguide,  multiple 
sizes  of  waveguide  and  corresponding  material  samples  must  be  used.  Using  a  different 
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geometry  of  waveguide,  such  as  Double  Ridge  Waveguide  (DRWG),  can  provide  a  larger 
bandwidth  of  wave  propagation  from  the  same  guide  [3,  11].  The  increased  bandwidth 
can  provide  higher  resolution  for  detection  of  defects  in  RAM,  thus  providing  a  more  cost 
effective  tool  used  in  coating  diagnosis.  An  example  of  both  geometries  is  shown  in  Figure 
1.1.  While  replacing  the  rectangular  waveguide  with  DRWG  will  expand  bandwidth,  new 
challenges  such  as  sample  fabrication  and  NWA  calibration  are  introduced.  Additionally, 
the  mathematical  modeling  of  a  DRWG  is  more  difficult.  To  ensure  accuracy,  an  in  depth 
analysis  of  a  NWA  calibration  scheme  will  need  to  be  investigated  as  well  as  exploring 
material  fabrication  challenges. 


Figure  1.1:  Double  Ridge  Waveguide  (left)  and  Rectangular  Waveguide  (right) 


1.2  Limitations  and  Challenges 

A  limitation  of  using  DRWG  for  materials  measurements  is  the  same  as  any  other 
waveguide  system;  the  bandwidth  will  still  be  limited  [15].  Even  though  a  single 
DRWG  can  significantly  increase  the  measurement  bandwidth  compared  to  rectangular 
waveguides,  measurements  could  be  desired  in  a  much  larger  band.  Additionally,  low 
frequency  measurements  are  extremely  difficult  to  attain  using  waveguide  methods,  as  the 
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waveguide  dimensions  become  very  large.  Two  challenges  arise  as  a  direct  result  of  the 
complex  shape:  sample  fabrication  and  mathematical  modeling  of  wave  propagation  in 
the  guide.  Samples  are  placed  in  a  holder  that  has  the  same  shape  and  dimensions  as  the 
waveguide  itself.  The  sample  should  fill  the  waveguide  completely,  so  the  geometry  for 
DRWG  is  clearly  more  complex  than  that  of  rectangular  waveguide  and  some  materials 
may  not  be  easily  formed  to  the  correct  shape.  The  mathematical  description  of  the 
electromagnetic  boundary  conditions  give  rise  to  a  much  more  intense  analysis,  necessary 
to  determine  the  constitutive  parameters. 

A  practical  limitation  of  taking  measurements  with  this  method  is  that  it  is 
fundamentally  destructive  in  nature.  While  the  sample  being  measured  may  not  be 
damaged,  the  shape  must  be  harvested  from  the  original  test  article.  For  example:  If  a 
measurement  of  an  aircraft’s  RAM  coating  is  desired,  the  sample  will  need  to  be  cut  from 
the  aircraft,  thus  damaging  the  LO  treatment.  This  work  will  provide  a  proof  of  concept  in 
order  to  further  research  non-destructive  techniques  that  are  extremely  useful  in  LO  aircraft 
diagnosis. 

1.3  Scope 

Waveguides,  both  rectangular  and  double  ridged,  can  be  purchased  in  a  variety 
of  different  sizes  and  dimensions.  In  inventory  at  AFIT  are  rectangular  waveguide 
measurement  devices  that  operate  in  the  X-Band  range  (8-12GHz).  To  accomplish  relevant 
comparisons  between  the  two  different  waveguides,  a  frequency  band  of  about  6-18GHz 
has  been  chosen  to  prove  the  concept.  This  frequency  range  is  commercially  available 
as  a  WRD-650  double  ridged  waveguide.  This  range  will  span  a  portion  of  the  C-band 
(5.85-8.2GHz)  and  include  all  of  the  X,  and  Ku  (12.4-18GHz)  bands  of  frequencies. 
Measurements  will  be  compared  among  rectangular  and  DRWG  for  validation,  even  though 
rectangular  waveguide  requires  three  bands  to  span  the  entire  DRWG  bandwidth. 
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1.4  Research  Goals 


As  stated,  the  primary  goal  of  this  research  is  to  acquire  the  ability  to  make  broadband 
material  measurements  from  a  single  material  sample.  The  broader  goal  is  to  save  time 
in  performing  these  measurements  and  reduce  the  amount  of  samples  necessary.  New 
measurement  tools,  such  as  MATLAB  code  to  externally  calibrate  data,  would  be  important 
to  future  work  in  this  field.  Investigation  on  material  shape  effects  in  DRWG  is  desired, 
since  the  existing  requirement  is  difficult  to  fabricate  compared  to  traditional  rectangular 
waveguide  samples.  The  uncertainty  and  effects  of  imperfect  material  samples  will  be 
investigated  as  well.  While  the  mathematical  model  of  DRWG  has  been  addressed  [15,  20], 
it  is  the  goal  of  this  thesis  to  perform  a  more  descriptive  method.  The  primary  output  of 
the  mathematical  methods  that  this  thesis  is  interested  in  is  the  cutoff  wavenumber  (kc)  or 
frequency  (/c),  which  are  related. 

1.5  Resource  Requirements 

Comparing  the  results  of  using  DRWG  to  rectangular  waveguide  will  verify  that  the 
methods  used  are  valid.  To  make  this  comparison,  it  will  be  necessary  to  have  all  the 
equipment  required  to  accomplish  both  types  of  measurements.  The  resources  necessary  to 
complete  the  objective  related  to  this  thesis  are  as  follows: 

•  MATLAB  software  for  external  calibration, 

•  Network  Analyzer  to  take  measurements, 

•  Rectangular  waveguide  with  associated  coaxial  cable  adapter  (2), 

•  Double  ridged  waveguide  with  associated  coaxial  cable  adapter  (2), 

•  Double  ridged  line  standard  to  verify  methodology, 

•  Coaxial  cabling  with  SMA  connectors  to  connect  the  measurement  guides  (2), 
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•  Metal  short  plate  for  ends  of  waveguides  for  reflect  calibration, 

•  Material  samples  in  rectangular  shape  for  comparison, 

•  Material  samples  in  double  ridged  shape  for  verification  of  methodology. 

1.6  Organization 

This  thesis  is  divided  into  five  chapters.  The  second  chapter  will  go  over  basic 
electromagnetic  and  material  concepts,  scattering  parameters,  waveguide  propagation 
fundamentals  and  will  include  an  overview  of  the  NWA  with  associated  calibration 
schemes.  The  third  chapter  will  derive  the  cutoff  frequency  for  DRWG,  explore  the 
Nicholson-Ross-Weir  (NRW)  method  for  determining  the  constitutive  parameters,  and 
discuss  sample  preparation.  The  fourth  chapter  will  discuss  results  relating  to  the 
constitutive  parameters,  uncertainty,  sample  modifications  and  comparisons  to  rectangular 
waveguide  measurements.  The  fifth  chapter  will  summarize  the  work  performed  for  this 
thesis  and  discuss  any  changes  or  improvements  for  future  work. 
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II.  Background 


To  better  appreciate  the  ability  to  take  measurements  of  the  constitutive  parameters  of  a 
material,  a  conceptual  knowledge  of  the  principles  involved  is  important.  Becoming 
familiar  with  existing  measurement  methods  is  important  to  understanding  the  methods 
presented  in  this  thesis.  This  chapter  will  discuss  electromagnetic  properties,  material 
concepts,  waveguide  fundamentals,  alternate  measurement  techniques  and  the  calibration 
method. 

2.1  Basic  Electromagnetic  Principles  and  Maxwell’s  Equations 

Electromagnetics  is  the  study  of  systems  of  charge  in  a  general  state  of  motion  [3]. 
While  many  scientists  have  contributed  to  the  modern  field  of  electromagnetic  theory,  the 
basis  equations  for  study  are  known  as  Maxwell’s  equations.  They  describe  the  behavior 
of  electromagnetic  waves  and  are  represented  in  either  integral  of  differential  form.  While 
examination  of  Maxwell’s  equations  is  broad  in  scope,  this  thesis  will  focus  on  specific 
attributes  and  can  be  summarized  succinctly.  In  the  following  sections,  important  concepts 
such  as  duality,  electric  and  magnetic  fields,  mono-poles  and  dipoles  will  be  reviewed. 
Some  common  symbols  used  in  this  thesis  are: 

•  E  -  Electric  Field  Intensity  (volts/meter) 

•  H  -  Magnetic  Field  Intensity  (amperes/meter) 

•  D  -  Electric  Flux  Density  (coulombs/square  meter) 

•  B  -  Magnetic  Flux  Density  (webers/square  meter) 

•  J  -  Electric  Current  Density  (amperes/square  meter) 

•  9Jt  -  Magnetic  Current  Density  (volts/square  meter) 
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These  symbols  are  denoted  as  vector  quantities  and  are  assumed  to  be  time  varying. 

2.1.1  Duality. 

In  mathematics,  if  two  equations  have  the  same  form,  their  solutions  will  also  have  the 
same  form.  This  concept  is  called  duality.  To  demonstrate  this  principle,  a  comparison  of 
the  Maxwell-Faraday  equation  and  Ampere’s  Law  can  be  used: 


v  x  e  =  -m  -  — 

dt 

(2.1) 

dD 

V  x  H  =  J  +  — . 

ot 

(2.2) 

Observing  that  (2.1)  and  (2.2)  have  similar  forms,  duality  can  be  applied.  This  principle 
is  important  in  the  study  of  electromagnetics  as  one  can  have  the  solution  to  one  type  of 
problem,  and  interchange  the  symbols  to  have  the  solution  to  another  problem  [11]. 

The  premise  of  using  duality  in  Maxwell’s  equations  assumes  an  existence  of  magnetic 
current  density  (9JI)  and  thus  magnetic  charges.  Physically  these  entities  have  not  been 
shown  to  exist,  but  can  be  of  great  assistance  in  the  analysis  and  study  of  electromagnetics 
[3].  Physically,  an  electric  charge  does  exist.  Gauss’s  Law  for  D  and  B  is  written  as, 


V  •  D  =  qev 

(2.3) 

(2.4) 

where  qev  is  electric  charge  and  qmv  is  the  magnetic  charge.  In  most  cases,  the  magnetic 
charge  can  be  set  to  zero  although  sometimes  it  is  beneficial  to  leave  it  as  a  variable  so 
duality  or  equivalence  principles  can  more  easily  be  observed  and  applied  [18]. 

2.2  Types  of  Materials 

Atomic  composition  determines  the  type  of  interaction  a  material  will  exhibit  in  the 
presence  of  an  electric  or  magnetic  field.  The  characteristics  of  this  interaction  describe 
the  properties  of  the  material.  These  properties  are  important  to  consider  when  analyzing 
a  phenomenon  or  attempting  to  design  a  product.  The  primary  categories  explored  in  this 
thesis  are  conductors,  dielectrics  and  magnetics,  but  others  exist  [3]. 
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2.2.1  Conductors. 


The  outer  shell  of  an  atom  is  known  as  the  valence  shell.  Electrons  that  occupy  this 
shell  are  known  as  valence  electrons.  For  electromagnetic  interaction  purposes,  this  shell 
is  of  prime  interest.  If  there  are  unmatched  valence  electrons,  they  are  known  as  “free” 
electrons.  If  a  material  is  comprised  of  molecules  with  a  lot  of  free  electrons,  it  is  called 
a  conductor.  In  conductors,  the  electrons  can  move  in  random  directions  and  different 
velocities  with  no  field  applied.  Movement  of  electrons  is  known  as  current  and  with  no 
field  applied,  the  net  current  will  be  zero  in  the  conductor  [3].  When  a  field  is  applied,  these 
free  electrons  will  move  with  a  direction  and  velocity  relative  and  proportional  to  the  field 
vector.  The  measure  of  how  susceptible  a  material  is  to  supporting  an  electric  or  magnetic 
current  is  called  Conductivity  and  is  denoted  by  cr. 

2.2.2  Dielectrics. 

As  opposed  to  conductors,  an  ideal  dielectric  does  not  have  free  electrons  available  to 
move.  Instead  of  conducting  a  current  in  the  presence  of  an  applied  field,  these  materials 
distort  or  change  shape.  This  resultant  distortion  is  known  as  a  dipole.  While  dipoles  do 
have  a  field  associated  with  them,  they  are  randomly  oriented  and  distributed  throughout 
a  material  and  will  create  a  net  field  of  zero.  When  a  field  is  applied  to  an  ideal  dielectric 
the  dipoles  may  align  themselves  with  the  field  vector,  known  as  polarization,  but  will  not 
allow  transfer  of  electrons.  The  measure  of  how  susceptible  a  material  is  to  becoming 
polarized  is  called  Permittivity  and  is  denoted  by  e  [1 1]. 

2.2.3  Magnetics. 

Magnetics  are  essentially  the  dual  of  dielectrics  mathematically.  Physically,  they  are 
different  for  the  same  reasons  described  in  Section  2.1.1.  Instead  of  a  simple  polarization 
and  corresponding  electric  field,  magnetic  materials  experience  a  torque.  Thinking  of 
dipoles  as  electric  current  loops  instead  of  dipoles  can  help  to  understand  this  concept 
as  shown  in  Figure  2.1.  The  torque  follows  the  right  hand  rule  similar  to  the  concept  of 
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moment  of  inertia  in  the  physics  of  motion.  When  the  moment  vectors  all  point  in  the  same 
direction,  the  material  is  said  to  be  magnetized.  The  measure  of  how  susceptible  a  material 
is  to  becoming  magnetized  is  called  Permeability  and  is  denoted  by  jj  [3]. 


&  c 

<±> 

j 

A) 

t  * 

CD 

Figure  2.1:  Magnetic  Loop  Illustration:  No  Magnetic  Field  (left)  and  Magnetic  Field 
Applied  (right) 


2.3  The  Constitutive  Parameters 

The  conductivity,  permittivity  and  permeability  of  a  material  are  collectively  known 
as  the  constitutive  parameters.  These  parameters  characterize  the  interaction  of  electric 
and  magnetic  fields  with  the  material.  When  applying  Maxwell’s  Equations,  (2.1)  through 
(2.4),  the  constitutive  parameters  can  be  used  to  make  some  relations.  The  equations  that 
demonstrate  the  relations  are  called  the  “Constitutive  Relations”  [11]  and  are  written  as, 


J  =  crE 

(2.5) 

D  =  eE 

(2.6) 

B  =  yUH. 

(2.7) 

Since  the  materials  studied  here  will  be  classified  largely  as  dialectics,  only  the 
permittivity  and  permeability  will  be  emphasized.  Further  use  of  the  term  “constitutive 
parameters”  will  imply  only  these  two  parameters.  Because  they  are  meant  to  describe 
the  medium  in  which  the  field  exists,  permittivity  and  permeability  can  be  decomposed  to 
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differentiate  between  free  space  and  a  material: 


e  =  ere0 

(2.8) 

H  =  AFjU  o, 

(2.9) 

where  er  and  n,  are  permittivity  and  permeability  relative  to  the  material  and  e0  and  fJo 

are  indicative  of  free  space.  The  free  space  parameters  have  the  approximate  and  defined 

values 

e0  =  8.854  x  1(T12  (farads/meter)  (2.10) 

juq  ~  4tt  x  10  7  (henry/meter).  (2.11) 

The  permeability  of  free  space  is  defined  to  be  the  above  value.  It  can  be  related  to  the  free 
space  permittivity  by 

6o  =  (2.12) 

c2H  o 

where  c  is  the  measured  speed  of  light  (~  2.99792  m/s). 

For  the  analysis  performed  in  this  work,  the  constitutive  parameters  are  not  directly 
measured.  Instead,  the  scattering  (or  S)  parameter  data  is  taken  from  the  NWA.  Later  on, 
the  NRW  algorithm  is  used  to  derive  the  constitutive  parameters  from  the  S-parameters. 

2.4  Scattering  Parameters 

The  overall  goal  of  this  research  effort  is  to  find  the  constitutive  parameters  of  a 
material.  Because  the  network  analyzer  does  not  measure  these  parameters  directly,  it  is 
vital  to  understand  what  measurements  are  actually  taken.  The  S-parameters  are  measured 
over  a  spectral  range  by  the  network  analyzer.  By  conservation  of  energy,  a  wave  can 
either  transmit  through,  be  absorbed  into,  or  reflect  from  a  material  [11].  A  graphical 
representation  of  reflection  and  transmission  is  shown  in  Figure  2.2.  A  simple  method 
of  describing  this  interface  is  by  utilizing  the  scattering  matrix.  It  is  possible  to  relate  the 
waves  transmitted  and  reflected  on  the  right  side  in  Figure  2.2,  to  the  waves  on  the  left 
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C2 

b2 


side  with.  The  relation  is  accomplished  with  the  scattering  parameters  [8].  A  mathematical 
representation  of  a  material  interface  is  shown  by 


by  -  C\S  11  +  b2S  12 
C2  =  CyS  21  +  bnS  22- 


(2.13) 

(2.14) 


Rewriting  in  a  matrix  form, 


by 

S  11 

S  12 

Cl 

C2 

S  21 

S  22 

b2 

(2.15) 


where  S  n  and  S  22  are  the  reflection  coefficients  from  the  right  and  left  sides  respectively. 
5 21  and  S 12  are  the  transmission  coefficients  from  the  right  and  left  sides  respectively. 
These  coefficients  are  collectively  referred  to  as  the  S-parameters,  and  are  represented  as 
complex  numbers  to  indicate  magnitude  and  phase. 

Observation  of  (2.15)  shows  that  the  left  side  of  the  matrix  equation  does  not 
correspond  to  the  left  side  of  Figure  2.2.  To  simplify  analysis,  the  T-parameters  can  be 
introduced.  T-parameters  can  mathematically  depict  the  Figure  2.2  by. 


(2.16) 


by  T 11  T 12  Z?2 

C\  T2\  T2  2  c2 

The  conversion  from  S  to  T  parameters  is  discussed  in  Section  2.8.  One  of  the  advantages 
to  using  the  T-parameters  is  the  simplicity  in  representing  a  network  consisting  of  multiple 
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layers.  The  T- matrices  can  be  written  as  matrix  multiplication  with  a  T- matrix  representing 
each  layer.  For  example,  if  Figure  2.2  included  n  number  of  blocks  representing  different 
mediums,  then  the  system  could  be  represented  mathematically  as: 

T  'T'T  rr  1  rr  1  ^2  ^2  'T'n  'T'n 

11  1  12  _  1  11  1  12  1  11  1  12  1  11  1  12 

T  rpj  y  1  rj-<2  T  U  TU 

21  1  22  J  L  21  i22jL  21  1  22  J  L  21  1  22 

2.5  Waves  and  Waveguide  Propagation 

To  better  understand  the  discussion,  a  brief  overview  of  the  type  of  waves  is  provided 
here.  The  first  concept  is  wave  modes  and  polarization.  The  modes  discussed  will  include 
Transverse  Electromagnetic  (TEM),  Transverse  Electric  (TE)  and  Transverse  Magnetic 
(TM).  A  TEM  wave  is  described  a  wave  where  both  E  and  H  are  contained  in  a  plane 
transverse  to  the  direction  of  propagation.  This  type  of  wave  can  exist  in  free  space  or 
in  a  guide  with  two  separate  conductors.  A  TE  wave  refers  to  the  E  field  in  the  plane 
transverse  to  the  direction  of  propagation.  Finally,  a  TM  wave  contains  the  H  field  in  the 
plane  transverse  to  the  direction  of  propagation  [3]. 

In  order  to  solve  for  the  constitutive  parameters,  the  cutoff  wavenumber  (or  frequency) 
is  needed  for  the  NRW  discussed  later.  A  cutoff  frequency  is  the  frequency  below  which 
a  wave  ceases  to  propagate  in  a  waveguide  for  a  given  propagation  mode  [11].  This  thesis 
is  exploring  DRWG  as  a  measurement  apparatus,  but  the  technique  is  similar  to  using 
rectangular  waveguides.  The  mathematical  goal  of  finding  the  cutoff  wavenumber  is  fairly 
straightforward  with  regard  to  rectangular  waveguide.  For  this  analysis  it  can  be  assumed 
that  the  wave  will  be  a  TE1  wave,  meaning  it  is  a  TE  wave  with  propagation  in  the  z- 
direction.  The  E  and  H  fields  are  desired,  but  are  extremely  difficult  to  solve  for  directly. 
It  is  common  practice  to  use  the  vector  potentials  A  (magnetic  vector  potential)  and  F 
(electric  vector  potential)  to  aid  in  the  analysis  [4].  For  this  analysis, 


ii 

ni> 

(2.18) 

A  =  0 

(2.19) 

(2.17) 


T 

T 
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where  Fz  is  the  z  component  of  the  potential.  For  any  propagating  wave  in  a  source  free 
region,  the  Helmholtz  Wave  Equation  must  be  satisfied 

V2FZ  +  JcF-  =  0  (2.20) 


where  k  is  the  wavenumber.  For  a  wave  standing  in  the  x  and  y  directions  and  traveling  in 
the  z  direction, 

Fz  =  (A  cos  kxx  +  B  sin  kxx){C  cos  kyy  +  D  sin  kyy)(Ee~jkzZ  +  FejkzZ)  (2.21) 


where  A.B.C.D.E,  and  F  are  constants  to  be  solved  for  and  j  =  V-I.  Now  that  the  vector 
potentials  are  known,  the  E  and  H  fields  can  be  solved  for.  The  equations  relating  the  fields 
to  the  vector  potentials  are 

1  T  1  V  x  F 

E  = -  k2 A  +  V(V  •  A) -  (2.22) 

jcoe/u  [  J  6 

and 

H  =  +  -1—  k2 F  +  V(V  •  F)  .  (2.23) 

H  jcoen  [ 

Unique  solutions  will  be  dependent  on  boundary  conditions  and  wave  mode.  The  boundary 
conditions  can  be  solved  for  in  terms  of  the  mode  numbers.  While  an  infinite  number  of 
modes  can  exist,  the  lowest  order  (or  dominant)  mode  for  rectangular  waveguide  is  known 
as  the  “1-0”  mode,  notated  as  TEzl0.  Detectors  in  the  NWA  used  in  this  work  detect  this 
dominate  mode,  so  it  is  of  particular  interest  [3]. 

2.5.1  Rectangular  Waveguide. 

Using  the  assumption  that  a  waveguide  is  constructed  using  a  Perfect  Electric 
Conductor  (PEC)  material  and  it  has  dimensions  a  and  b.  the  cutoff  frequency  can  be  shown 


to  be  given  by 


(2.24) 


where  n  and  m  describe  the  mode  [11].  A  convenient  property  of  rectangular  waveguide  is 


the  ability  to  express  the  cutoff  frequency  in  closed  form.  As  derived  in  later  analysis,  this 
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is  not  possible  for  DRWG.  In  addition  to  cutoff  frequency  information,  it  is  important  to 
understand  how  the  fields  are  distributed  in  the  cross  section  of  the  guide. 

2.5.2  Field  Behavior. 

For  rectangular  waveguide,  the  electric  field  inside  the  guide  can  be  described  as 
sinusoidally  distributed  along  the  x-axis,  oriented  in  the  y-direction  [3].  With  respect  to 
material  measurements  performed  in  rectangular  waveguide,  the  sample  should  fill  the 
entire  guide  completely. 

For  DRWG,  the  fields  are  highly  concentrated  at  the  corners  of  the  ridge  [20].  While 
measurement  error  will  occur  if  the  material  does  not  completely  fill  the  guide,  it  stands  to 
reason  that  large  errors  would  be  expected  if  any  gaps  exist  in  the  vicinity  of  the  corners  of 
the  ridge.  This  error  will  be  explored  in  a  later  chapter. 

2.6  Measurement  Bandwidth 

As  stated,  the  primary  advantage  to  using  DRWG  over  rectangular  waveguide  is  the 
increased  bandwidth.  Table  2.6  shows  an  example  of  manufacturer  supplied  bandwidth 
capabilities  for  some  example  waveguides  [14]  [6].  As  can  be  seen,  it  will  take 
approximately  three  rectangular  waveguides  to  accomplish  the  same  measurement  as  one 
DRWG.  Some  verification  of  this  information  will  be  explored  in  later  chapters. 


Table  2.1:  Waveguide  Bandwidth 


Waveguide  Type 

Model 

Bandwidth 

Rectangular  (C-Band) 

WR-137 

5.85-8.20 

Rectangular  (X-Band) 

WR-90 

8.20-12.40 

Rectangular  (Ku-Band) 

WR-62 

12.40-18.00 

DRWG 

WRD-650 

6.5-18.00 
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2.7  Other  Measurement  Methods 


This  thesis  focuses  on  waveguide  measurements  that  support  TE  or  TM  propagation 
modes  only.  There  are  more  techniques  that  can  be  utilized  to  derive  the  constitutive 
parameters  while  performing  similar  measurement  methods,  with  respect  to  similar 
algorithms.  Three  other  methods  include  stripline,  focus  beam  and  coaxial  techniques. 
The  devices  used  in  these  alternate  methods  can  support  a  TEM  wave  in  addition  to  TE  or 
TM  modes. 

2.7.1  Coaxial  Cable. 

Coaxial  cables  are  used  commonly  in  a  vast  array  of  applications  due  to  their  ease  of 
use,  variety  of  connectors,  relatively  robust  design  and  most  importantly,  large  bandwidth. 
This  type  of  transmission  line  consists  of  a  center  conductor  surrounded  by  a  dielectric 
material  and  then  encased  in  another  conductor  [5].  An  example  is  shown  in  Figure  2.3. 

Because  of  the  ability  to  transmit  TEM  waves,  the  bandwidth  of  coaxial  cabling  can  be 
extremely  large.  The  lower  limit  is  theoretically  DC  while  the  upper  limit  is  determined 
by  the  excitement  of  higher  order  modes.  One  major  drawback  to  using  coaxial  cable  for 
constitutive  parameter  measurements  is  the  discontinuity  of  the  radial  electric  field  [2].  If 
an  air  gap  exists  between  the  center  conductor  and  the  sample  under  test,  the  uncertainty  in 


Outer  Conductor 


Inner  Conductor 


Figure  2.3:  Diagram  of  Coaxial  Cable 
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the  measurement  grows  large.  Coaxial  cable  measurements  can  also  suffer  from  inaccuracy 
due  to  poor  field  penetration  into  the  material  [5]. 

2.7.2  Stripline. 

Stripline  is  a  type  of  planar  transmission  line,  analogous  to  coaxial  cable  without  the 
outer  conductor  completely  surrounding  the  dielectric  material,  as  seen  in  Figure  2.4.  It  is 
constructed  of  two  plates  with  a  conductor  running  lengthwise  in  the  center.  Also  similar 
to  coaxial  cable,  stripline  supports  TEM  wave  propagation  modes  and  is  band-limited  for 
the  same  reasons.  The  method  used  to  make  material  measurements  in  a  stripline  is  to  form 
a  material  sample  around  the  center  conductor.  Again  there  is  a  significant  challenge  to 
designing  a  material  sample  to  fit  this  area.  Gaps  between  the  conductors  (inner  or  outer) 
and  the  material  will  cause  uncertainty  in  measurement  similar  to  coaxial  cable. 


Sample _ 

\ 

/ 

Inner  Conductor 

Outer^Conductor 

(a)  Diagram 


(b)  Picture 

Figure  2.4:  Stripline  Diagram  and  Picture 
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2.7.3  Focus  Beam. 


The  focus  beam  system  consists  of  hom  antennas,  lenses  and  a  sample  holder 
to  make  measurements  (Figure  2.5).  The  horn  antennas  launch  a  spherical  wave  to 
the  lenses,  that  in  turn,  transform  the  wave  into  a  plane  wave  [10].  Because  of  the 
equiphase  and  equiamplitude  characteristics,  plane  waves  make  deriving  the  intrinsic  wave 
properties  easier,  thus  producing  a  measurement  from  the  material  interaction  [3,  10,  16]. 
Sample  fabrication  is  more  easily  accomplished  for  the  focus  beam  system  with  the 
trade  off  being  complications  of  taking  a  measurement.  Some  additional  advantages  of 
the  focus  beam  system  is  the  ability  to  take  off  normal  or  oblique  measurements,  and 
the  ability  to  heat  the  sample  since  it  doesn’t  touch  any  waveguide  [5].  Calibration  of 
the  system  requires  the  lenses,  the  sample  and  sample  holder  to  exist  at  exact  known 
locations.  Additionally,  scattering  from  the  sample  holder  and  other  “clutter”  sources  can 
be  introduced  into  the  measurement  because  the  system  is  not  contained  like  coaxial  cable, 
stripline  or  waveguide  [10].  Accurate  measurements  require  extremely  good  control  of  the 
environment  surrounding  the  system  and  accurate  calibration  techniques.  Sample  sizes  can 
also  be  a  limiting  factor.  The  bandwidth  is  determined  by  the  physical  dimensions  of  th 
lenses  (lower  limit)  and  the  antennas  or  cabling  (upper  limit). 


Figure  2.5:  Picture  of  Focus  Beam  Measurement  System 
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2.8  Network  Analyzer  Operation  and  Calibration 

A  rough  diagram  of  a  materials  measurement  configuration  is  shown  in  Figure 
2.6(a).  The  detectors  that  perform  the  actual  measurements  are  contained  inside  the  NWA. 
Everything  that  electrically  connects  the  detectors  to  each  other  must  be  mathematically 
accounted  for.  A  block  diagram  characterizing  this  system  is  shown  in  Figure  2.6(b). 


Cal  1  Cal  1 


(a)  Network  Analyzer 


Det  1  Call  Call  Cal  2  Cal  2  Det  2 


(b)  NWA  Block  Diagram 

Figure  2.6:  Network  Analyzer  Test  Configuration 
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Note  that  Figure  2.6(a)  shows  “Port  1”  and  “Port  2”.  These  ports  are  not  the  detectors, 
only  where  the  coaxial  cable  connects  to  the  device.  Calibration  will  adjust  for  all  of  the 
components  up  to  the  calibration  planes,  notated  as  “Cal  1  and  Cal  2”  in  Figure  2.6(b).  The 
electrical  path  from  a  detector  to  a  calibration  plane  (Figure  2.6(b))  is  known  as  the  error 
adapter  [19].  There  are  two  error  adapters  that  will  need  to  be  represented,  “A”  and  “B”  for 
reference  here.  The  block  “S”  is  the  sample. 

The  system  can  be  represented  by  the  matrix  equation, 


(2.25) 


where  the  superscript  ms  is  the  detector  measured  value,  A,  B  and  S  are  represent  the  blocks 
from  Figure  2.6(b).  The  T  matrix  can  be  related  to  the  S -parameters  by 


1 


T  =  — 
52i 


|  S 21 S  12  -  S  nS  22  s  11  | 

-S  22  1 

Rewriting  (2.25)  in  terms  of  (2.26)  and  using  superscripts  to  describe  each  block, 
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1 

cA 
o  21 


cA  cA  _  cA  cA  oA 
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CB 
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'YB 


Rearranging  and  identifying  the  unknowns, 


(2.26) 


(2.27) 


(2.28) 


^A  'j1# 

it  is  shown  that  there  are  seven  unknown  values  to  solve  for:  S  S  ,  S  2I S  An,  S  21 5  f2,  S  Au , 
Sfj,  S*2,  and  SBr 

2.8.1  TRL  Calibration. 

Now  that  the  seven  unknowns  are  identified,  solutions  can  be  derived.  It  is  necessary 
to  perform  three  calibration  measurements  to  find  solutions.  The  first  is  the  “Thru” 
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(Figure  2.7(a))  measurement  which  is  comprised  of  only  the  two  waveguides  being  coupled 
together.  The  second  is  the  “Reflect”  measurement  which  involves  placing  a  solid  metal 
block  on  the  end  of  each  guide,  where  the  calibration  planes  are  located.  The  reflect 
measurement  is  shown  in  Figure  2.7(b)  with  two  reflection  standards  placed  between  the 
waveguides,  thus  electrically  terminating  each  guide.  The  final  calibration  measurement  is 
the  “Line”  (Figure  2.7(c))  which  requires  the  placement  of  a  line  standard  (also  known  as 
the  sample  holder)  in  between  the  waveguides.  In  addition  to  recording  the  S -parameter 
values  for  each  of  these  measurements,  the  “port  mis-match”  terms  will  also  be  collected 
from  the  NWA.  These  terms  quantify  the  amount  of  power  collected  from  each  port 
reflecting  from  the  opposite  port  and  are  displayed  as  a  ratio. 


(a)  Thru 


(b)  Reflect 


(c)  Line 

Figure  2.7:  TRL  Calibration  Standard  Configurations 
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2.8. 1.1  Thru  Measurement. 


Writing  (2.28)  using  superscripts  to  indicate  the  “Thru”  measurement  and  utilizing  the 
short  hand  indicated. 


1 


c  A  c  B 
°21°21 


(2.29) 


results  in  a  matrix  equation.  Here,  the  focus  is  the  T'  term.  For  a  perfect  waveguide,  there 
should  be  no  reflection  and  all  the  energy  will  be  transmitted.  Written  mathematically, 


S[l,St22  =  0  (2.30) 

5^,5i2=l.  (2.31) 


Therefore,  the  value  of  the  Tr  term  will  be  the  identity  matrix.  This  allows  for  solutions  to 
the  terms  contained  in  the  matrix, 


T'  = 


o  21 


O  t  Of  PI  Of 

°  21°  12  °  22 

cf 
^  11 
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_  cf 

O  22 

1 

0 

1 

(2.32) 


The  two  terms  with  the  tilde  over  them  are  unknowns  that  can  be  solved  for  using  the  Thru- 
Reflect-Line  (TRL)  measurements,  after  the  common  coefficient  is  factored  out.  Solving 
(2.29)  for  Tb  yields, 

Tb  =  TA~'TmtS2lS2i-  (2-33) 


Now  that  an  equation  has  been  derived  for  one  of  the  common  tilde  denoted  terms,  the 
derivation  will  proceed  to  the  “Line”  measurement. 

2.8. 1.2  Line  Measurement. 

In  addition  to  physically  holding  the  test  sample,  the  line  standard  also  serves  to  have 
a  standard  phase  measurement.  The  standard  itself  is  simply  a  short  section  of  waveguide 
that  aligns  well  to  the  rest  of  the  system.  Applying  the  same  notation  as  in  section  2.8. 1.1 
for  the  line  measurement. 


rpm/  _ 


c  A  c  B 
°21°21 


(2.34) 
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(a)  Rectangular  Standards:  Line  (L)  and  Re-  (b)  DRWG  Standards:  Line  (L)  and  Reflect 
fleet  (R)  (R) 

Figure  2.8:  Waveguide  Standards 


For  a  perfectly-aligned  line  standard,  the  S-parameters  associated  with  reflection  will  be 
zero  and  the  transmission  parameters  will  simply  undergo  a  phase  change.  Mathematically, 
this  can  be  written  as 


r./  n/ 

o  1 1 ,  a  22 

W  rW 


V'  S1  -  p-jkzl  -  p 

°  21’  °  12  —  ^  —  r  ? 


0  (2.35) 

(2.36) 

where  the  l  in  the  exponential  is  the  length  of  the  line  standard,  P  is  used  for  simplicity  when 
re-writing,  and  kz  is  wavenumber  associated  with  the  propagating  wave  in  the  z-direction 
(down  the  guide).  Following  the  same  procedure  as  section  2. 8. 1.1,  T1  can  be  written  as 

1 


T'  = 


Substituting  (2.33)  into  (2.34)  yields, 


p 2 

0 

P  0 

0 

1 

o 

(2.37) 


rjV7?/  _  rji/irji/rjvl  1  rjv/if 


(2.38) 


By  rearranging  the  terms  in  (2.38),  another  variable  can  be  introduced  to  aid  in  solving  the 
matrix  equation.  This  is  shown  by 


rjVi 


(2.39) 
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where  M  is  a  matrix  representing  the  product  of  Tm/  and  T"“  1 .  Rearranging  terms  again  to 
eliminate  the  inverse  matrices  produces 


MTa  = 

(2.40) 

Writing  (2.40)  in  matrix  form  and  substituting  in  (2.37)  produces 

mn  mn  ^12  ^12 

m2 1  m22  An  A22  An  A-22 

P  0 

(2.41) 

0  P"1 

where  the  m  and  A  entries  are  elements  of  M  and  TA  matrices  respectively.  By  performing 

matrix  multiplication  on  (2.41)  produces  four  equations: 

m\  1  An  +  m\2A2\  -  P An 

(2.42) 

ninAn  +  m22A2i  =  PA2 1 

(2.43) 

mn  An  +  m\2A22  =  P  ^  An 

(2.44) 

mnAu  +  mnAu  =  P  *A22. 

(2.45) 

Since  there  is  a  common  term  (P)  in  the  first  two  equations,  the  next  step  will  be  to  write 
one  of  them  begin  equivalent  to  P  so  it  can  be  substituted  in.  Rearranging  (2.43)  results  in, 

A\ 

m2  i  — — h  m22  =  P.  (2.46) 

A2i 

Substitution  of  eq  (2.46)  into  eq  (2.42) 

_  _  f  An 

mnAu  +mnA2i  =  m2l—  +  m2l  An,  (2.47) 

v  ^21  , 

yields  an  equation  that  resembles  a  quadratic.  Dividing  (2.47)  through  by  A2i  and 
rearranging  provides  a  quadratic  equation: 

'A  'f2  A 

mu  +  (m22  -  mn)  ^  -  ml2  =  0.  (2.48) 

vA2i7  A2i 
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Following  the  same  process,  but  for  (2.44)  and  (2.45),  produces  a  second  quadratic 


equation: 


m2 1 


'A12' 


J 


+  (m22  -  mu)  -  mi2  =  0. 

A22 


(2.49) 


Equations  (2.48)  and  (2.49)  both  satisfy  the  same  form  using  the  quadratic  formula, 
An  An  _  -  (m22  -  mn)  +  V(m22  -  mn)2  +  4m2iffli2 

A21  A22  2/W21 

Applying  the  matrix  entries  that  constitute  An,  A2i,  An,  and  A22, 

T  cA  cA  _  cA  c/ 

An  J  21°  12  °  11 J  ; 


(2.50) 


11 
A21 


-SA 

0  22 


A  cAl 

22  _  ~A 

-  0  n 


c  A  c  A 
O  2|  O  12 


cA 

o  22 


=  (2 


(in  _  -2, 

a22-  1 


(2.51) 


(2.52) 


where  a  and  b  are  used  for  simplicity  in  rewriting.  For  a  typical  network  analyzer,  the 
reflections  should  be  low  since  there  is  just  an  empty  guide,  so  the  terms  in  (2.51)  can  be 
related  using  an  inequality, 


|sfi| 


« 


r  A  r  A 

o  21o  12 


cA 

o  22 


(2.53) 


Using  this  information,  it  is  determined  that  b  is  the  small  root  and  a  is  the  large  root  of  eq 
(2.50).  By  subtracting  b  from  a  and  simplifying. 


b  -  a  =  S  j  j  -  S  fj  - 


cA  cA 
°  21°  12 

cA 
O  22 


cA  cA 
°  21°  12 
cA 
°  22 


(2.54) 


and  rearranging  for  one  of  the  seven  unknowns  mentioned  previously, 


sisn  =  5 22  (b  -  a) . 


(2.55) 


At  this  point,  expressions  have  been  found  for  two  of  the  seven  unknowns.  Now  it  will 
be  necessary  to  revisit  the  thru  and  line  measurement  expressions  again,  but  for  different 
terms. 
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2.8. 1.3  Recalculations  of  Thru  and  Line  Expressions. 

Using  (2.29)  and  (2.32)  and  rearranging  to  solve  in  terms  of  TB  will  provide  more 
equations.  By  making  some  observations  about  the  variables,  more  independent  equations 
can  be  realized.  Solving  for  TA  will  yield 

T4  =  Tm,fB-'S^Srf.  (2.56) 


Following  the  same  process  shown  above  (section  2. 8. 1.2),  a  matrix  expression  can  be 
written  as 


(2.57) 


N 

where  N  is  shorthand  for  the  product  Tmf  '  T",/.  Rearranging  and  writing  (2.57)  to  eliminate 
the  inverse  matrix  produces, 


Tbn  =  T  iTb 


(2.58) 


Since  individual  expressions  are  desired,  the  next  step  is  to  rewrite  in  matrix  terms  so 
multiplication  can  be  used; 


B  ii  B 12 
B2]  B2  2 


n  n  n  12 
n2\  n22 


P  0 
0  p-1 


B\\  B\2 

B2\  b22 


(2.59) 


Again,  explicitly  writing  out  the  matrix  multiplication  operation  will  produce  four 
equations: 


Bnnn  +  B\2n2\  =  PBn  (2.60) 

B\\U\2  +  B[2n22  =  PBn  (2.61) 

B2\n\\  +  B22n2\  —  P  * B2\  (2.62) 

B2iti\2  +  B22n22  =  P  ^ B22.  (2.63) 


Solving  (2.61)  for  P, 


Bn 

~  n\2  +  n22 

Bn 


=  P, 


(2.64) 
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and  substitution  into  (2.60)  gives, 


Bnitn  +  -  — — i-  ^22  \B{\. 

U12  ) 


(2.65) 


Rearranging  and  dividing  by  B\2  results  in, 


Bn)  ,  B 11 

—  ni2  +  (n22  ~  nn)  — - n2l  =  0. 

Bn)  Bn 


(2.66) 


Using  the  same  procedure,  but  for  (2.62)  and  (2.63),  another  equation  of  similar  form  as 
(2.66)  is  formed, 


(£21)  ^  ,b21 

—  nn  +  (n22  -  nn)  — - n2i  =  0. 

V  B22  /  b22 

Finally,  a  quadratic  equation  similar  to  (2.50)  is  realized, 

Bu_  Bn.  _  ~  (n21  -  nn)  +  V(n22  -  nn )2  +  4n2Rh2 
B\\  B22  2h  12 

Showing  the  two  solutions  as  products  of  their  constitutive  S -parameters, 


(2.67) 


(2.68) 


c  B  c  B  c  B 

0  12  °  11°  22  __ 

c  B  c  B 

CB  °  21°  12 

(2.69) 

cB  ~ 

O  22 

^22  ~C 

3  11 

B21  _  ~S22  _ 
B22  1 

-S  22  =  d. 

(2.70) 

where  two  more  variables,  c  and  d,  are  used  for  simplicity  in  rewriting.  Again,  using  the 
fact  that  a  typical  NWA  will  have  low  reflection, 


c  B  c  B 

I  cB  |  ^  °21°  12 
|  '  22 1  «  „B 
O  22 


(2.71) 


So,  c  will  be  the  large  root  and  d  will  be  the  small  root  of  (2.68).  Subtracting  d  from  c  and 
rearranging  will  yield  a  solution  for  another  of  the  seven  original  unknowns 


S^SBn  =  Snu  (c-d). 


(2.72) 


Four  of  the  original  seven  unknowns  have  been  found.  Additional  measurement 
expressions  will  be  required  to  find  the  remaining  terms. 
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2.8. 1.4  Reflect  Terms. 

The  reflect  standard  can  be  shown  mathematically  from  Port  1  as  [19], 

cA  cA  rr 

c mr  _  cA  ,  °21°  121 

S”=5ii+1-S£F- 


(2.73) 


Solving  for  P  yields 


(c mr  _  cA  \ 

°11 


(2.74) 


(cAcA  cApA\’ 

where  the  superscript  mr  signifies  the  measured  reflect  terms.  Substituting  in  the  variables 
already  defined,  a  and  b.  (2.74)  can  be  written  as, 


S'"[-b 


(2.75) 


s  ihlT  -4 

For  a  short  that  is  placed  with  no  offset  from  the  end  of  the  waveguide,  the  value  of  T' 
should  be  extremely  close  to  -1,  indicative  of  a  strong  reflection.  Applying  similar  steps 
starting  at  (2.73)  for  the  reflect  parameter  from  Port  2, 


n  mr 
^  22 


oi  rfi  rr 
r.B  J  12J  211 
0  22  + 


i-SfiP 


and  substituting  in  c  and  d. 


rr  = 


S%  +  d 


S  f,(ss  +  c)' 

will  result  in  a  similar  equation  to  (2.77).  Setting  (2.75)  equal  to  (2.77), 


(2.76) 


(2.77) 


S%  +  d 


S™-b 

sUsu~a)~  su  bs+9 


(2.78) 


and  solving  for  S  22  yields 


r*  A  _  r»  B 

°  22  “  °  11 


65+96n  -*) 


(2.79) 


6s+46'.t  -«)’ 

Five  of  the  seven  unknowns  now  have  expressions.  For  the  final  two,  a  revisit  of  the  thru 
expressions  will  have  to  occur.  The  following  expressions  assume  that  there  is  no  detector 
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mis-match.  This  is  not  true  for  a  real  NWA,  so  later  on  correction  for  the  mis-match  will 
be  derived.  If  the  detector  mis-match  has  been  corrected,  then 


amt  _  a  A 

°  11  “  °  11 


a  A  a  A  a  B 
°21°  12°  11 

1  _  rA  rB  ■ 
r  o  22o  j  j 


(2.80) 


Solving  for  S  fj  and  substituting  in  a  and  b  gives, 


cB  _ 

o  ii  - 


sn;\-b 


rA 
O  22  | 


which  is  the  sixth  unknown.  Applying  (2.81)  to  (2.79)  provides  an  expression  for  S^: 

(ST!  -  b)  (i?2r  +  c)(sn  ~b) 


(2.81) 


CA  _ 

O  22  - 


’  22 


(s";l-a)(s^  +  d){sz-a] 


(2.82) 


Multiplying  both  sides  of  the  equation  by  S  22  and  taking  the  square  root  of  the  entire 


expression  results  in 


-  + 

°22  _  - 


(s’«-i>)(s^-b)(sT  +  c) 


fl/2 


(2.83) 


(s~-a)(S7[-a)(s%  +  d)_ 

and  a  root  choice  is  required.  Using  (2.75)  and  the  fact  that,  for  a  reflect  standard 
F  w  -1+0 j,  the  real  part  of  F  should  be  less  than  zero.  To  use  this  information,  each  root 
of  S 29  will  have  to  be  substituted  into  (2.75)  and  the  appropriate  root  chosen.  To  find  the 
seventh  and  final  unknown,  the  S  2|  will  be  revisited  to  provide  an  additional  expression. 


Solving 


c*  mt 
O  2\ 


oA  c A 
o2i02i 


1  _  VA  cB  ’ 
1  022On 

and  rearranging  for  the  final  unknown  shows, 


(2.84) 


a  A  a  B  _  r*  mt  / 1  _  a  A  a  B  \ 

°21°21  “  °21  \L  °22°  U)  ‘ 


(2.85) 


Now  all  of  the  unknowns  have  been  found,  but  the  port  mis-match  problem  has  not  yet  been 
addressed. 
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2.8. 1.5  Port  Mis-match  Correction. 

From  Rytting  [19],  the  source  and  load  mismatch  terms  can  be  determined  from  the 
following  matrix  equation: 


h  W  cm  cm  n 

lm  Im  _  11  12  hn  lm 

W  C m  C m  r.  n' 

u2m  u2m  ^  21  u  22  u2m  Li2m 


(2.86) 


Rearranging  to  solve  for  5'"  terms, 


C  III  C  I 
O  U  J, 


cm  c  in 

P  21  J  22  J 


Tm  °\m\ 


' 2m  u2m 


K  a\ma2 m  aima2m  _a2m 


(2.87) 


Performing  matrix  multiplication  and  simplifying. 


c  m  _ 
O  n  - 


erf  _  c  d  aim 

J11  ^  12  aim 


aim  a^m 


c  III  _ 
0  12  _ 


c  m  _ 
0  22  _ 


a  I  m  ^2m 

erf  _  erf  aim 

J21  J  22  flj, „ 

J  aim  ^  lm 

a  lm  alfn 

erf  _  e  rf  17 lm 

J22  J21a' 


6^1/n 


(2.88) 


(2.89) 


(2.90) 


(2.91) 


where  the  superscript  d  indicates  the  detector  measurement.  The  a  and  b  terms  are 
available  as  additional  output  from  the  NWA,  displayed  as  ratios,  R2/RI  and  RI/R2.  The 
measurements  correspond  to  the  equations  above  by, 


R2  eh... 

Port  1,  —  =  — 

/v  1  @\ m 


Port  2, 


R}_  =  <n 

R2  a2m 


(2.92) 


(2.93) 


All  the  terms  necessary  to  correct  for  port  mis-match  have  been  found  and  the  TRL 
calibration  scheme  is  complete. 
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III.  Methodology 


The  previous  chapter  demonstrated  the  algorithms  and  mathematics  required  to 
understand  a  NWA  based  measurement  system.  The  equipment  configuration 
using  DRWG  and  data  acquired  will  appear  similar  to  those  of  rectangular  waveguide. 
Calibration  will  be  accomplished  in  MATLAB  using  the  four  measurements  detailed  below. 
The  final  step  is  to  derive  the  constitutive  parameters  using  the  NRW  algorithm  [17,  21]. 

3.1  Double  Ridged  Waveguide  Analysis 

To  make  the  analysis  of  the  DRWG  tractable,  the  geometry  can  be  broken  down 
into  three  regions  and  the  common  coordinate  system  will  be  located  in  the  center  of  the 
guide,  as  shown  in  Figure  3.1  [12].  The  point  at  the  upper  right  corner  of  region  I  will  be 
considered  (Ax,  Ay).  The  point  located  at  the  upper  right  of  region  II  is  referred  to  as  (|,  |J. 

The  first  step  for  this  analysis  will  be  to  define  the  fields  using  the  vector  potentials,  F 
and  A.  Since  the  coordinate  system  has  been  chosen  to  have  the  z  direction  as  the  waveguide 
axis  and  an  electric  source  is  assumed  (TEZ), 


ii 

(3.1) 

A  =  0. 

(3.2) 

The  forcing  function  for  the  propagating  wave  in  the  z  direction  is  given  by 

Fz  =  (B  cos  kxx  +  C  sin  kxx)  ( D  cos  kyy  +  E  sin  kyy)  (A+e~jkzZ  +  A~e'kzZ)  (3.3) 

fix)  g(y)  h+(z ) 

where  A.B.C.D  and  E  are  constants  and  k  is  the  wavenumber  in  a  given  direction.  For 
simplicity  and  conciseness,  the  terms  f(x),  g(y )  and  h+(z)  will  be  used  in  place  of  the 
terms  in  (3.3).  The  electric  field  intensity  (E)  and  the  magnetic  field  intensity  (H)  can  be 
expressed  as 

E  =  — —  \k2 A  +  V(V  •  A)1  -  (3.4) 

jcoe/j.  1  1  6 
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H  =  +  — —  \k2 F  +  V(V  •  F)1 . 

H  juen  L  J 

Applying  (3.1)  and  (3.72)  into  (3.4)  and  (3.5)  yields 

V  x  F 


E  =  — 


H  =  — \ic F  +  V(V  •  F)1 . 
u,)(=n  l  J 


J0J6/U 


(3.5) 

(3.6) 

(3.7) 


To  continue,  each  field  vector  will  need  to  be  broken  up  into  its  respective  field  components. 
This  is  accomplished  by  performing  the  operations  in  (3.6)  and  (3.7).  The  field  components 


are: 


Ex  =  -j  f(x)(D  sin  kyy  -  E  cos  kyy)h+(z) 

(3.8) 

Ey  =  -—(Bcoskxx  -  C  sin  kxx)g(y)h+(z) 
e 

(3.9) 

Ez  =  0 

(3.10) 

k  k 

Hx  -  (B  sin  kxx  -  C  cos  kxx)g(y)  (A+e~jkzZ  -  A~ejkzZ) 

h~(z) 

(3.11) 

k,k~ 

Hy  -  ——  f(x)(  D  sin  kvy  -  E  cos  kxy)h  (z) 
jcoen 

(3.12) 

Hz  =  —(k2  +  k2z)Fz 
jcoen  > - - - - 

(3.13) 

=  .  f(x)g{y)h+{z) 

jcoen 

(3.14) 

where  kc  is  the  cutoff  wavenumber  of  the  guide.  Now  that  the  general  fields  have  been 
derived,  analysis  can  continue  by  deriving  the  terms  for  each  region.  Since  the  boundary 
conditions  for  each  region  will  vary  from  each  other,  each  region  will  be  derived  separately 
and  the  coordinate  system  will  have  an  origin  at  the  center  of  the  guide.  The  configurations 
for  analysis  are  shown  in  Figure  3.1. 


3.1.1  Region  1. 

Region  I  contains  a  PEC  at  the  top  and  bottom  boundaries  which  are  located  at  +Ay. 
This  region  is  also  called  the  “gap”  region.  Since  the  origin  of  the  coordinate  system  is 
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(a)  Coordinate  System  (b)  Regions 

Figure  3.1:  DRWG  Configurations  for  Analysis 


chosen  to  be  at  the  center  of  the  guide,  there  will  be  a  “shift”  term  substituted  in  for  the 
variables  in  each  term.  For  region  I,  the  field  equations  shown  generically  in  (3.8)-(3.14) 
can  be  written  as: 


(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 

(3.21) 


where  (y  -  Ay)  has  been  substituted  in  for  x.  At  the  upper  and  lower  PEC  boundaries 
(y  =  ±Ay),  the  tangential  electric  fields  will  be  zero.  In  this  case,  the  x  and  z  directed  fields 
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will  be  zero  and  can  be  written  as, 


Exx  =  0.  (3.22) 

y=±Ay 

2=(-oo,+oo) 

Applying  the  boundary  condition  to  (3.15)  will  yield  two  equations: 

EXl(y  =  Ay)  =  D{  sin  kyi (0)  -  E\  cos  kyi(0)  =  0  (3.23) 

EXl(y  =  -Ay)  =  Di  sin  kyt (-2Ay)  =  0.  (3.24) 

Solving  the  two  equations  will  yield  a  result  for  one  of  the  constants  and  for  the  y  associated 
wavenumber: 


E  i=0 
nn 

kyi  =  2Ay ' 


(3.25) 

(3.26) 


Since  there  is  no  PEC  boundary  on  the  right  or  left  side,  more  analysis  will  need  to  be  done 
before  applying  a  boundary  condition  to  the  fields  there. 

3.1.2  Region  II. 

This  section  will  analyze  region  II  as  defined  in  Figure  3.1(b).  This  region  will  also 
be  referred  to  as  the  “trough”  region  in  later  notation.  Again,  the  assumption  that  the 
waveguide  walls  can  be  modeled  as  a  PEC  has  been  made.  The  forcing  function  will  be  of 
the  same  form  as  in  the  previous  analysis  and  will  have  a  different  “shift”  term  associated 
with  it.  The  function  for  region  II  is  given  by, 

h+(z ) 

riOO  g2(y ) 

(3.27) 


Z2 


LI  LI  (/  I 

B2  cos  kX2  -  -  x  +  C2  sin  kXl  -  -  x  D2  cos  ky2  y  -  -  +  E2  sin  k. 


W2 


y_2 
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where  h+(z )  is  the  same  as  before  and  the  subscript  “2”  signifies  region  II.  The  field 
equations  are  given  by: 

EX1  =  ky,f2(x)  [d2  sin  kyi  |y  -  ^ j  -  E2  cos  kyi  Py  -  ^)J j  h+(z)  (3.28) 

Ey 2  =  ~kX2  {l B2  sin  kX2  _  *)  “  C2  cos  kX2  ^  -  x)  jj  g2(y)h+(z )  (3.29) 

EZ2  =  0  (3.30) 

Hx ,  =  [Bi  sin kX2  -  C2  cos  kX2  ^  -  ^)j)g2Cv)/^(z)  (3-31) 

H}2  =  ico/u^2^  (D2  Sin  kyi  r  _  _  El- cos  k*2  _  ^)J j  h~(z)  (3.32) 

k2 

HZ2  =  Mx)g2(y)h+(z )  (3.33) 

J^H 

where  h~(z)  is  the  same  as  in  section  3.1.1.  The  PEC  exists  at  the  x  value  and  y  values 
±|,  which  can  be  written  as 


=  0 

(3.34) 

V- 

II 

1+ 

MO 

£=(— 00,+00) 

=  0. 

(3.35) 

y=i 

£=(— 00,+00) 

Applying  the  PEC  boundary  conditions  to  the  field  equations  results  in 


£2  =  0  (3.36) 

C2  =  0  (3.37) 


and 


(3.38) 

(3.39) 

By  inspection  of  Figure  3.1(b),  region  II  and  region  III  will  be  similar  in  analysis. 


D2  sin  ky2(b)  =  0 
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3.1.3  Region  III. 

The  geometry  of  regions  II  and  III  will  be  similar,  but  the  variables  will  differ  due  to 
different  locations  of  the  PEC  boundary.  Region  III  is  also  called  a  “trough”  region.  Again, 
starting  with  the  forcing  function,  modified  for  region  III, 

Fz 3  =  (fi3  COS  kX3  ^  +  *1  +  C3  sin  kX3  ^  +  .vjj  (d3  cos  ky3  |y  -  ^  j  +  £3  sin  ky3  (y  -  ^ \ \  h+(z ) 


gi(y) 


The  field  equations  are  given  by: 


-  KJ-iix)  ( D3  sin  ky3  ( y  -  -  I  -  £3  cos  ky3 1  y  -  -)  | )  h+(z) 


Ey3  =  ~kx 3  ( B3  sin  kX3  +  v)  -  C3  cos  /cA;i  +  w) ) )  g3(y)h+(z) 


£7,  =0 


31  X3  = 


Hy3  = 


^3  sin  Cr3  (^+  X)~C3  cos  kX3  +  *)))  #3  00^  (z) 

(D3  sin  (y  -  -  £3  COS  C3  (>’  -  ^)Jj  /C(z) 


= 


-f3(x)g3(y)h'(z) 


J^n 

The  boundaries  are  located  at  y  =  ±|  and  jc  =  -f .  Shown  mathematically, 


(3.40) 

(3.41) 

(3.42) 

(3.43) 

(3.44) 

(3.45) 

(3.46) 


-,v3 


0 


v-+t 

y-±2 

Z=(-o o,+oo) 

Ey3  =  0. 

y=-i 

Z=(-oo,+oo) 

Applying  the  PEC  boundary  conditions  to  the  field  equations  results  in 


(3.47) 

(3.48) 


£3  =0 

C3  =0 


(3.49) 

(3.50) 
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and 


D3  sin  ky3(b)  =  0 


(3.51) 


v3 


m7r 


(3.52) 


Additionally,  comparing  the  expressions  for  the  wave  numbers  in  regions  II  and  III,  the 
observation 

kX2  =  kx  3  (3.53) 


can  be  made.  Equation  (3.53)  will  be  useful  in  determining  the  rest  of  the  unknown 
equations  for  this  analysis. 

3.1.4  Continuous  Fields  Boundary  Conditions. 

The  final  boundary  condition  to  apply  is  the  existence  of  continuous  fields  at  the 
transitional  line  between  each  region.  This  region  will  be  referred  to  as  the  “gap”  boundary 
region.  The  boundary  conditions  require  that  the  tangential  electric  and  magnetic  fields  be 
continuous  at  this  gap  [15].  For  this  analysis,  let 


Ei(y')  be  the  gap  electric  field  at  x  =  A.v  (3.54) 

E2(y')  be  the  gap  electric  field  at  x  =  -Ax.  (3.55) 


Writing  the  field  expression  and  substituting  in  the  values  found  in  Sections  3.1.1  -  3.1.3, 
the  expressions  for  E’i(y')  and  E2(y')  can  be  written  as  [15], 


E\{y)  =  ~Y  kx  1„  ( BUt  sin  kXlnx  -  Cln  cos  kXlnx)  cos  j  ^  (y'  -  Ay) J  = 


Z,  .  ,  a  \  mn  / 

kX2mB2m  sm  kX2m  -  -  A*  I  cos  —  y 


(3.56) 


E2(y')  =  ^  sin  kXlnx  +  Cu  cos  kXlnx)  cos  j  YY  (/  -  by) 


Y  kXlm^m  sin  kX2m 


mn 


Ax  cos  —  y' 


(3.57) 
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A  series  representation  of  the  fields  will  be  difficult  to  solve  for  the  unknowns.  Performing 
integration  on  the  expression  will  “filter”  the  series  to  acquire  the  “nth”  term.  From  the 
identity  [13], 


a+L 


I 

a-L 


mnx  nnx 
cos  — — -  cos  -—ax 


0 
L 
2  L 


m  =£  n 
m  =  n  ±  0 
m  =  n  =  0 


(3.58) 


it  can  be  observed  that  the  series  can  be  simplified  by  multiplication  the  original  function 
(3.56)  and  integration  over  the  boundary.  Application  of  (3.58)  yields  the  following  two 
equations: 


Av 


/ 

-Ay 


E\(y')  cos  iy' 
2Ay 


Ay)dy'  =  ~kXhi(B 


'In 


sin  kx,  Ax 

-*1  n 


C\n  cos  kx lwAx)  A y  (1  +  ^,o) 


(3.59) 

/YY17T  I  b\  I  Cl  \  b 

E\(y  )  cos  —  1/  -  -Jdy'  =  kX2mB2m  sin  kX2m  -  Ax j- (l  +  Smfi) ,  (3.60) 

-Ay 


where  the  Kronecker  Delta  function  facilitates  the  solution  given  in  (3.58).  Similar 
expressions  involving  E2(y')  can  be  written, 


Ay 

/TUT 

E2(y')  cos  —  ( y '  -  Ay)  dy  =  kXh 

-Ay 


(B i„  sin  kXh  Ax  +  Cln  cos  kXll  Ax)  Ay  (1  +  6tu0) 


(3.61) 

Ay 

/WITT  I  b\  l  a  \  b 

E2{y)  cos  —  Iy'  -  -  J  dy'  =  -kX2mB3m  sin  kX2m  -  A.vj  -  (1  +  6mfi) .  (3.62) 

-Ay 
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Solving  for  B2m  and  B3m  in  (3.60)  and  (3.62)produces 

Ay 

f  W)cosf  (y'-l)dy' 

- Ay 

~  — h - : - 7 - \ 

kX2m2^  ^m,o)  ^ i n  (2  AxJ 

Ay 

-  f  e2(/) cos  f  (/-I) dy 

-Ay 

B’im  =  - - - 7 - T- 

kX2m  f  (1  +  c>m,o)  sm^  (f  -  A.i'j 


(3.63) 

(3.64) 

(3.65) 


To  simplify  further  analysis,  the  numerators  in  (3.63)  and  (3.64)  will  be  referred  to  as  I3 
and  I4  respectively.  If  (3.59)  and  (3.61)  are  rearranged  to 


Ay 

f  E1(y)cos^(y-l)dy 

—Ay 

Cu,  cos  kXlnAx  -  Bu  sin  kX]nAx  =  — - - — — — - — — - 

kxin^y  \*-  On, 0/ 

Ay 

f  E2(y)  cos  %(y-l)dy 

—Ay 

c in  cos  kXl  Ax  +  B  1„  sin  kXlnAx  =  — - - — —— - — — - , 

kXliAy(l  +6nfi) 


(3.66) 

(3.67) 

(3.68) 


it  will  be  possible  to  solve  for  C ln  and  BUl  with  simple  expressions  if  the  same  shorthand 
is  used  as  before.  The  numerators  of  (3.66)  and  (3.67)  will  be  referred  to  as  Ii  and  I2 
respectively.  Now  C\„  and  BUl  can  be  written  as 


_ Ii  +  h _ 

2kXlnAy(l  +  Snfi)  cos  kXhAx 

_ h^h _ 

2 kXln  Ay  ( 1  +  Snfi)  sin  kXln  Ax  ‘ 


(3.69) 

(3.70) 


Now  that  these  constants  have  been  found,  (3.56)  and  (3.57)  can  be  written  by  substituting 
in  the  values  found  in  (3.69)  and  (3.70).  The  gap  magnetic  fields  must  by  continuous  by 
the  boundary  condition  [15],  so  equating  the  tangential  magnetic  field  at  the  boundary  will 
produce  continuity.  For  this  analysis,  the  z-directed  fields  will  be  used.  Selecting  the  field 
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in  region  I  and  equating  it  to  the  field  in  region  II  at  the  gap, 


Hz  \{x  =  -Ax)  =  Hz  2(x  =  -Ax), 


(3.71) 


the  following  expression  is  formed  after  substitution  of  the  shorthand  and  constants 
previously  solved  for: 

V  7i  (tan  kXln  Ax  -  cot  kXln  Ax)  +  I2  (tan  kXln  Ax  +  cot  kXln  Ax)  tm 


2kXhAy{\  +<S„,0)) 


ms—(y.Xy) 


,  I, cot kx2m  t,  A.vJ  mn  l  b\ 

/  - t - cos  y  -  -  =  0. 

^  k,^(l+Sm,0)  hV  2 1 


(3.72) 


Applying  the  same  technique  and  boundary  condition  to  the  gap  field  between  region  I  and 


region  III, 


Hz  i(x  =  -Ax)  =  Hz  $(x  =  -Ax), 


(3.73) 


a  similar  expression  as  (3.72)  will  arise: 

Zl i  (tan  kXln  Ax  +  cot  kXln  Ax)  +  I2  (tan  kXlnAx  -  cot  kXln  Ax)  nn 

2kXlnAy  (1  +  6nfi))  2  Ay 

„  I 4  cot  kX2m  -  A.vj  mn  /  /A 

V  kX2J_(l+6my)  b  X  2) 


(3.74) 


Since  the  goal  for  this  analysis  is  to  obtain  the  cutoff  frequency  for  use  in  later  calculations, 
solving  for  the  wavenumbers  will  be  critical. 

3.1.5  Finding  the  Cutoff  Frequency. 

As  stated,  the  primary  purpose  of  this  analysis  for  this  work  is  to  find  the  cutoff 
wavenumber.  For  the  DRWG  geometry,  the  cutoff  wavenumber  can  be  related  to  the  other 
wavenumbers  associated  with  each  region  by, 


]P-  -  lr 2  +  h2 
K-c  ~  X-yg  -I-  Xxg 


]p-  -  Jr  2  +  ]A 

Kc  Kyt  ^  Kxt 


(3.75) 

(3.76) 


where  kc  is  the  cutoff  wavenumber,  and  the  subscripts  g  and  t  indicate  the  gap  and  trough 


modes  [15].  The  cutoff  wavenumber  for  this  system  is  the  point  when  kz  is  zero,  so 
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generation  of  equations  that  will  result  in  zero  for  kz  is  important.  To  do  this,  E\(y'),  E2(y'), 
Ii,  I2,  I3  and  I4  will  be  rearranged  and  the  unknown  modal  amplitudes  (a  and  J3)  will  be 
substituted  in, 


Ei(y  )  =  Xa/C0S(2Av)^ 

(3.77) 

Eiiy  )  =  Yj cos  (  2Av )  ^  A-y 

(3.78) 

h  =  ^  o'/ Ay  ( 1  +61,0)  6  i,n 

1 

(3.79) 

h  =  ^  /7/ Ay  ( 1  +6i,o)6[,n 

1 

(3.80) 

h  =  ^  «,iA+(m,  /) 

/ 

(3.81) 

/4  =  2 PiV(m,l ) 

(3.82) 

1 


where 


>P+(  m,  l ) 


1 

2 


k  +  k 

*vln  T  Kyi n 


sin  |  mn  \  - 


(—  1  )n  sin  \mn\  i  ^ 


(3.83) 


To  continue  the  analysis,  the  testing  function  is 


p(x,  y)dy 


(3.84) 


where  p(x,y )  are  the  field  expressions  (eqs  (3.72)  and  (3.74))  derived  previously. 
Substitution  of  eqs  (3.77)-(3.82)  into  the  field  expressions  and  the  testing  function,  the 
following  equations  are  derived: 


tan  kxU,Ax  -  cot  kxln Ax 
eta - zv - (1  +  6a, 0)  Ay+ 


2k 


xlh 


tankxmAx  +  cotkxUAx 
Pn - - (1  +  6a,  0)  Ay- 


27: 


xlh 


EE 


if/+(m ,  ri)ifj+{m,  h)  cot  kx2m  (f  -  A*) 

k-x2m  2  (1  +  6m_  0) 


0  (3.85) 
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tan  A*  +  cotk^Ax 

a,-, - — - (1  +  4,0 )  Ay+ 


— ^  \  i  // 


tannin  Ax  -  cot  kx\hAx 
fir - tt - (1  +  4,o)  Av- 


2k 


xlh 


„  „  (A  (m,  n)ijj  (nu  n)  cot  kx2m  (  f  -  Ax) 

>  >  /3n - t - — - -  =  0.  (3.86) 

nm  kx2m\  (1  +  Sm,o) 

Finally  a  matrix  equation  can  be  formed  from  eqs  (3.85)  and  (3.86)  and  can  be  written  as, 


A(kc)  [a i  ...  an  J3i  ...  fin]T  =  0 


(3.87) 


where  A(kc)  is  a  matrix  of  basis  columns  and  testing  rows  and  0  is  specifically  the  zero 
vector.  The  final  step  is  to  find  the  values  of  kc  that  satisfy  detA  =  0.  The  method  used  is 
to  find  the  eigenvector  of  A(kc),  which  will  satisfy  the  equation.  For  this  thesis,  Newton’s 
Method  was  used  by  guessing  at  an  initial  value  and  iteration  until  the  determinate  that 
was  equal  to  zero  was  found.  MATLAB  was  utilized  to  numerically  solve  for  the  cutoff 
wavenumber  or  the  dominant  and  next  order  mode  to  establish  bandwidth. 


3.2  Measurement  Procedure 

The  measurement  procedures  for  DRWG  and  rectangular  waveguide  are  identical 
in  practice.  While  order  is  not  particularly  important,  the  “thru”,  “line”,  and  “reflect” 
measurements  described  in  Section  2.8.1  need  to  be  taken  to  calibrate  the  system. 
Additionally,  the  sample  under  test  needs  to  be  measured  by  the  system.  The  physical 
thicknesses  of  the  line  standard  and  sample  also  need  to  be  recorded  for  use  in  the  NRW 
algorithm  described  next,  although  it  is  possible  to  determine  the  line  thickness  with  the 
thru  and  line  measurements.  Upon  collection  of  the  data,  processing  will  be  accomplished 
in  MATLAB. 

One  area  that  has  not  been  addressed  is  the  selection  of  the  line  standard.  It  is  desired 
to  have  a  large  sample  holder  to  maximize  possible  sample  thickness,  but  the  width  of  the 
holder  will  be  have  to  be  small  enough  to  not  cause  a  phase  ambiguity.  The  first  step  is 
to  design  the  length  so  that  it  is  approximately  one  quarter  of  the  wavelength  at  the  center 
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frequency.  The  next  step  is  to  verify  that  the  phase  does  become  ambiguous  in  the  line 
standard  when  a  wave  is  propagating  through  it  [9].  Figure  3.2  shows  a  7  mm  and  9  mm 
line  standard  phase  plot.  As  can  be  seen,  the  phase  for  7  mm  does  not  reach  a  singularity 
across  the  bandwidth  as  the  9  mm  standard  does.  To  maximize  sample  thickness  and  to 
keep  a  safe  margin  from  the  zero  phase  point,  the  7  mm  standard  was  chosen. 


Frequency  (GHz) 


Figure  3.2:  Phase  Plot  Depicting  7  mm  and  9  mm  Line  Standards 


3.3  Nicholson-Ross-Weir 

The  mathematical  technique  used  to  derive  expressions  for  the  constitutive  parameters 
is  called  the  NRW  algorithm.  This  technique  will  determine  the  constitutive  parameters 
from  the  S-parameters  relative  to  one  port  at  a  time.  For  homogeneous  and  isotropic 
media,  the  parameters  should  be  the  same  for  the  opposite  port  as  well.  For  the 
following  derivation,  the  S-parameters  involved  are  the  forward  parameters  and  represent 
the  calibrated  data,  denoted  as  Sexxxp  [5]. 
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The  first  quantity  desired  is  known  as  the  inter-facial  reflection  coefficient  ( R )  and  is 


given  by 


r  =  q±  ve2- 1, 


(3.88) 


where  Q  is  given  by 


(S“p)2  -  ( Se2\p )2  +  1 

rs  n  exp  ’ 

n 


(3.89) 


Equation  3.89  has  two  solutions  that  are  related  as  inverses.  One  solution  will  have  a 
magnitude  less  than  one  and  is  the  solution  used  to  proceed.  The  next  step  is  to  determine 
the  one-way  phase  delay,  denoted  as  P,  and  is  given  by 


1  -RSe* 


(3.90) 


The  cutoff  frequency  for  rectangular  waveguide  was  given  in  the  previous  chapter  (2.24). 
Applying  (2.24)  to  a  TE10  mode  wave,  the  cut  off  frequency,  k0  and  y0  are  given  by 


■j  ,  ko  =  ureoPo  .  To  =  K  ~  ko 


(3.91) 


where  a  is  the  long  dimension  of  the  waveguide  and  to  is  the  angular  frequency  (to  =  2  nf). 
The  final  step  in  the  NRW  algorithm  to  find  the  constitutive  parameters.  The  relative 


permittivity  ( er )  is  given  by 


k2c  -  ((In  P)/l)2 


b2]nP  l+R 
o  ro i  \  i -R 


(3.92) 


where  1  is  the  thickness  of  the  sample  [10].  The  relative  permeability  (jur)  is  then  computed 


c  In  P  / 1  +  R 


7o  l  U  -R 


(3.93) 


The  constitutive  parameters  can  be  determined  separately  for  forward  (5n,52i)  or  reverse 
parameters  (S22A 12),  by  simply  substituting  the  correct  reflection/transmission  parameters 
above  (S  22  for  S 1 1  and  S 12  for  S2i).  Since  the  calibration  planes  are  set  at  either  side  of 
the  sample  holder,  it  is  also  necessary  to  calculate  the  phase  delay  from  either  port  to  the 
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4 


(a)  S-parameters  calibrated  to  cal  planes  (b)  S-paremters  shifted  to  sample 

Figure  3.3:  S-Parameter  Calibration  and  Adjustment 


sample  and  back.  To  calculate  the  phase  delay,  the  path  to  and  from  the  sample  will  be 
accounted  for.  In  Figure  3.3(a),  the  calibration  planes  are  the  edges  of  the  holder.  The 
scattering  parameters  measured  are  shown  and  an  adjustment  must  be  made  to  account  for 
displacement  of  the  sample.  Starting  with  S ,  ]  and  S21, 

5,,  =e-Ao45-^e-Ao4  (3.94) 

S21  =  (3.95) 

where  kz0  is  the  propagation  wavenumber  in  free  space  and  the  other  variables  correspond 
to  the  values  in  Figure  3.3.  The  goal  is  to  obtain  the  S-parameters  for  the  sample,  so  solving 

r  rt  sample  j  sample 

ior  ij  |  j  an  q  13  21  9 

sample  _ 

sample  _  $  ^jkt o(4“4) 


(3.96) 

(3.97) 
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and  the  forward  parameters  are  found.  Following  a  similar  process  for  the  reverse 
parameters  will  show, 

Ss“2mpte  =  S22ej2k^  (3.98) 

Ss™ple  =  W*zo(4_4).  (3.99) 

The  phase  adjustment  for  the  S-parameters  has  been  adjusted  to  the  sample  itself  instead 
of  the  sample  holder  and  the  constitutive  parameters  can  be  extracted  using  the  NRW 
algorithm. 

3.4  Sample  Preparation 

For  this  work,  there  were  two  materials  measured.  The  first  was  a  plastic  material 
provided  by  Air  Force  Research  Laboratory  (AFRL),  an  acrylic  sample  and  two  different 
microwave  absorbers.  The  sample  constructed  by  AFRL  is  a  polymer  known  commercially 
as  VeroP  and  was  fabricated  in  a  Polyjet  3D  printer.  The  design  was  accomplished  in 
Solidworks  and  is  shown  in  Figure  3.4.  The  dimensions  used  to  draw  the  sample  were 
taken  from  datasheets  for  WRD-650.  The  dimensions  of  the  final  product  were  slightly 
bigger  than  the  guide.  Since  tight  samples  are  desired  for  maximum  accuracy,  the  samples 
were  forced  into  the  line  standard.  This  was  only  attempted  because  the  size  of  the  sample 
was  extremely  close  to  the  size  of  the  guide.  The  result  was  that  some  of  the  excess  material 
was  shaved  off  of  the  sample  by  the  line  standard  and  the  sample  fit  extremely  tight  in  the 
waveguide.  This  process  was  repeated  for  each  sample  with  similar  results.  In  order  to 
verify  the  designs  prior  to  sending  out  for  fabrication  by  AFRL,  an  acrylic  sample  was 
made  using  local  3D  printer  resources.  The  acrylic  was  able  to  fit  similarly  to  the  AFRL 
sample. 

The  absorber  material  was  formed  using  laser  cut  plastic  as  a  template,  then  cut  out 
using  a  razor.  The  same  procedure  of  pushing  the  line  standard  onto  the  sample  was 
performed.  The  original  thought  process  was  similar  to  that  of  the  more  rigid  plastic 
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Figure  3.4:  Sample  Drawing  (Values  in  Inches) 
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samples;  the  material  would  nicely  shave  off  and  fit  well  into  the  guide.  This  was  not 
true  in  practice,  especially  with  the  thicker  absorber  material  tested.  The  results  of  each 
method  will  be  discussed  and  testing  in  depth  in  the  next  chapter. 
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IV.  Results 


This  chapter  will  focus  on  the  extraction  of  the  constitutive  parameters  from  material 
samples.  All  of  the  measurements  were  performed  using  the  Agilent  E8362B 
NWA.  The  measurements  were  taken  in  the  frequency  domain  from  6-18GHz.  For 
each  measurement,  the  four  scattering  parameters  were  measured,  in  addition  to  the  port 
mis-match  power  ratios  at  each  detector.  External  calibration  was  conducted  using  a 
“Thru”,  “Line”  and  “Reflect”  measurement.  The  cutoff  frequency  is  also  determined 
using  MATLAB  code  and  the  bandwidth  of  DRWG  is  demonstrated.  Uncertainty  in  the 
measurements  will  be  presented  in  the  propagation  of  error  method.  Comparisons  to 
rectangular  waveguide  measurements  will  be  demonstrated  and  challenges  that  presented 
in  the  experimentation  will  be  discussed. 

4.1  Bandwidth 

From  Section  2.6,  a  table  of  usable  bandwidth  existed  for  DRWG  and  rectangular 
waveguide.  Table  4.1  shows  the  bandwidth  calculated  using  the  MATLAB  code  to  find  the 
cutoff  frequency. 


Table  4.1:  Calculated  DRWG  Bandwidth 


WRD-650 

Dominate  Mode  Cutoff 

Next  Order  Cutoff 

Bandwidth 

5.303  GHz 

18.291  GHz 

12.988  GHz 

Mathematically,  there  should  be  approximately  13  GHz  of  bandwidth  available  for  use 
from  DRWG.  In  reality,  a  safety  factor  needs  to  be  established  from  the  cutoff  frequency 
for  usable  bandwidth.  Figure  4.1  shows  a  vertical  line  at  the  calculated  cutoff  of  5.303 
GHz.  The  manufacturers  suggested  lower  band  is  also  shown  at  6.5  GHz.  While  6.5  GHz 
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Frequency  (GHz) 


Figure  4.1:  Cutoff  Frequency  Illustration  Using  Uncalibrated  Data 


is  free  from  errors  exhibited  by  the  lower  cutoff,  6.3  GHz  is  an  approximate  1.2x  safety 
factor.  Therefore,  the  measurement  bandwidth  presented  will  range  from  6.3  GHz  to  18 
GHz. 

4.2  Uncertainty  in  Measurements 

While  waveguides  can  provide  good  measurements,  uncertainty  will  manifest  in 
multiple  places.  The  technique  used  to  quantify  uncertainty  will  be  a  propagation  of 
error  technique,  since  a  statistically  significant  number  of  measurements  is  not  practical 
in  the  measurement  technique  used  here  [7].  The  types  of  uncertainty  that  will  arise  from 
the  system,  the  sample,  the  sample  position,  the  holder  and  the  waveguide  tolerances  are 
discussed.  This  section  will  also  display  results  with  uncertainty  factored  in  as  error  bars. 
All  of  the  uncertainty  values  used  are  displayed  in  Table  4.2.  Additionally,  all  of  the  values 
shown  represent  a  worst-case  scenario  for  uncertainty. 
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Table  4.2:  Uncertainty  Values 


Uncertainty  Type 

AFRL  Sample 

Acrylic 

SF-3.0 

FGM-125 

ISnl 

0.01125 

0.01125 

0.02 

0.02 

ZS  ii 

1.6° 

1.6° 

1.6° 

1.6° 

IS21I 

0.1  dB 

0.1  dB 

0.06dB 

0.06dB 

ZS  21 

0 

CO 

d 

O 

0 

0.5° 

0 

OO 

d 

Position 

0.01  mm 

0.01  mm 

0.05  mm 

0.10  mm 

Thickness 

0.01  mm 

0.01  mm 

0.05  mm 

0.05  mm 

Holder  Width 

0.035  mm 

0.035  mm 

0.035  mm 

0.035  mm 

kc 

5% 

5% 

5% 

5% 

4.2.1  System  Uncertainty. 

The  NWA  has  separate  published  uncertainties  for  reflection  and  transmission 
coefficients.  For  the  NRW  algorithm,  there  is  interest  in  the  phase  and  magnitude 
uncertainty  for  both  coefficients.  The  values  of  the  uncertainty  can  be  extracted  from  the 
users  manual  provided  on  the  NWA  as  plots  [1],  The  worst  case  values  were  used  for  the 
results  in  this  work.  In  Table  4.2,  the  system  uncertainty  is  the  first  four  rows  where  |5  n| 
and  ZSn  are  the  reflection  uncertainties  in  magnitude  and  phase  respectively.  |S2il  and 
ZS  2 1  denote  the  transmission  uncertainties  in  magnitude  and  phase. 

4.2.2  Sample  and  Sample  Holder  Uncertainty. 

The  sample  under  test  needs  to  be  measured  and  those  values  are  used  in  the  NRW 
algorithm.  Sample  thickness  and  dimensions  will  have  some  sort  of  uncertainty  associated 
with  them.  For  this  analysis,  only  the  sample  thickness  will  be  taken  into  consideration 
since  the  waveguide  will  determine  the  other  dimensions.  The  sample  thickness  varied 
for  the  type  of  material.  The  rigid  materials  (AFRL  sample  and  acrylic)  did  not  vary 
in  thickness  and  their  overall  uncertainty  was  set  at  the  lowest  possible  resolution  of  the 
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Figure  4.2:  A  Single  Sample  with  Two  Thicknesses 


calipers  used  to  make  the  measurement.  The  microwave  absorber  materials  were  not  as 
straightforward  to  measure.  In  some  cases,  a  single  sample  presented  a  gradient  across  the 
dimensions  as  seen  in  Figure  4.2. 

In  addition  to  the  sample  itself,  there  is  also  some  uncertainty  in  the  placement  of  the 
sample  in  the  waveguide.  When  installing  the  sample  into  the  holder,  efforts  are  made  to 
make  sure  the  face  of  the  sample  is  flush  with  the  face  of  the  guide.  The  sample  holder  is 
placed  on  a  flat  surface  and  the  sample  under  test  is  forced  down  so  the  displacement  from 
one  side  is  very  close  to  zero.  Again,  there  was  a  difference  in  uncertainty  for  the  ridged 
materials  and  the  softer  materials.  The  rigid  materials  presented  a  measured  uncertainty  of 
less  than  the  machine  precision  of  the  calipers,  while  the  softer  materials  varied  in  position 
slightly  more. 

The  sample  holder  itself  also  has  some  uncertainty  associated  with  it.  To  find  an 
uncertainty  term,  the  line  and  thru  measurements  were  used.  The  line  length  can  be 
calculated  by  using  the  phase  delay  found  by  comparing  the  two  measurements.  Recall, 


Ta  'mTa  = 


e~iht  o 
0  ejk*e 


(4.1) 
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from  Section  2. 8. 1.2.  By  finding  the  eigenvector  of  the  left  side,  an  overdetermined  solution 
will  result  and  €,  the  line  length,  can  be  solved  for.  This  results  in  an  electric  measurement 
over  the  bandwidth  of  interest,  which  is  shown  in  Figure  4.3.  Taking  the  standard  deviation 
of  the  data  plotted  is  the  uncertainty  value  used  here. 


Figure  4.3:  Line  Length  Calculated  Using  Thru  and  Line  Measurements 


4.2.3  Waveguide  Tolerance  Uncertainty. 

The  waveguide  itself  has  an  uncertainty  due  to  imperfect  manufacturing.  The  primary 
effect  of  imperfect  dimensions  on  uncertainty  is  in  the  cutoff  wavenumber,  kc.  While 
the  ideal  DRWG  will  have  completely  square  corners  on  the  ridges,  the  actual  physical 
specimens  are  not  square  due  to  limitations  in  machining  tools  and  fabrication  techniques. 
The  perfect  specimen  is  shown  in  Figure  4.4(a)  and  the  actual  product  is  shown  in  Figure 
4.4(b).  The  implications  of  the  imperfect  waveguide  will  impact  the  actual  value  of  kc, 
as  well  as  make  sample  fabrication  more  challenging.  The  value  by  which  kc  will  differ 
is  approximately  5%,  as  determined  by  another  technique  of  determination  [12].  The 
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Figure  4.4:  DRWG  Configurations  for  Analysis 


rounded  corners  will  actually  drive  the  cutoff  frequency  up,  which  will  also  have  the  effect 
of  decreasing  the  bandwidth  slightly. 

4.2.4  Application  of  Uncertainty  Terms. 

The  uncertainty  method  used  here  will  act  on  the  constitutive  parameter  equations. 
Each  uncertainty  term  will  be  calculated  and  added  in  quadrature.  The  formula  used  to 
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calculate  the  propagation  of  uncertainty  will  be. 


Reflection  Terms 


T ransmission  T erms 


^ - s  ^ ^ 


Magnitude  T erm  Phase  Term  Magnitude  T erm  Phase  Term 


Sample  Length  Line  Length  Cutoff  Position  Term 


(4.2) 


where  Sn  and  S2i  have  magnitude  and  phase  components  and  £,  d  and  k,  are  the  sample 
length,  line  length  and  cutoff  wavenumber  uncertainties  respectfully.  Similar  equations  are 
also  applied  for  the  imaginary  part  of  the  permittivity  (ej.)  and  both  parts  of  the  permeability 
(ju'r  and  n[). 


4.3  Measured  Sample  Results 

For  the  AFRL  sample,  the  measurements  with  uncertainty  are  shown  in  Figures  4.5(a) 
and  4.5(b).  As  can  be  seen  in  Figure  4.5(b),  the  relative  permeability  uncertainty  (jur) 
is  fairly  sensitive  to  changes  in  kc  at  the  lower  frequencies.  The  results  show  a  uniform 
response  across  the  bandwidth.  The  overall  uncertainty  is  low  across  the  bands  due  to 
small  sample  uncertainties  as  discussed  about  the  rigid  samples. 

During  the  sample  fabrication  process,  the  3D  printer  from  Air  Force  Institute  of 
Technology  (AFIT)  was  used  to  verify  the  sample  size.  The  material  used  is  a  commercially 
available  acrylic.  The  constitutive  parameters  for  this  material  are  shown  in  Figures  4.6(a) 
and  4.6(b).  The  results  show  that  the  material  has  a  fairly  uniform  response  across  the 
bandwidth  that  closely  resembled  the  sample  from  AFRL.  Again,  the  overall  uncertainty  is 
low  for  low  sample  uncertainties. 

Both  of  the  rigid  samples  above  have  very  similar  electrical  properties.  They  are  both 
polymers  that  were  made  on  similar  3D  printer  machines.  As  stated,  fabrication  was  made 
for  this  work  specifically.  Not  all  measurement  situations  will  allow  a  sample  to  be  easily 
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(a)  AFRL  Sample  Results  with  Uncertainty  ( er ) 


(b)  AFRL  Sample  Results  with  Uncertainty  (/ur) 
Figure  4.5:  AFRL  Sample  Results 


fabricated  by  precision  machinery,  so  attempting  to  measure  non-rigid,  custom  cut  pieces 
seems  to  be  a  more  realistic  application. 


54 


(a)  Acrylic  Sample  Results  with  Uncertainty  (er) 


(b)  Acrylic  Sample  Results  with  Uncertainty  (fjr) 

Figure  4.6:  Acrylic  Results 


Upon  initial  observations,  the  microwave  absorber  was  simple  to  cut  and  shape  into 
a  quality  measurement  specimen.  While  measuring,  it  was  quickly  discovered  that  the 
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(a)  SF-3.0  Absorber  Results  with  Uncertainty  (er) 


(b)  SF-3.0  Absorber  Results  with  Uncertainty  (/ur) 
Figure  4.7:  SF-3.0  Results 


fabrication  process  benefits  drastically  by  having  a  close  attention  to  detail.  The  easiest  of 
the  two  materials  used  was  the  SF-3.0  from  ECCOSORB.  It  is  thin  and  easy  to  cut,  making 
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(b)  FGM-125  Absorber  Results  with  Uncertainty  (//,  ) 
Figure  4.8:  FGM-125  Results 


it  straightforward  to  fabricate.  The  results  of  measurement  are  displayed  in  Figures  4.7(a) 
and  4.7(b).  The  final  material  measured  was  ECCOSORB  FGM-125,  a  commercially 
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available  microwave  absorber.  The  results  from  the  FGM-125  are  shown  in  Figures  4.8(a) 
and  4.8(b).  The  variance  and  difficulties  of  shaping  this  material  by  hand  are  discussed  in 
the  next  section. 

4.4  Sample  Modification 

Great  difficulty  was  encountered  in  fabricating  flexible  material  samples.  The  FGM- 
125  proved  to  be  the  most  difficult.  The  measurements  would  show  problems  in  the  higher 
frequency  range,  particularly  around  16  GHz  and  up.  Initially,  this  issue  was  falsely 
assumed  to  a  problem  discussed  in  Section  4.6  that  was  discovered  previously.  Upon 
observation  of  some  poorly  made  sample  post  technique  improvement,  it  was  discovered 
that  the  sample  is  very  sensitive  to  flaws.  Figure  4.9  shows  a  small  gap  on  one  side  of  the 
ridge  and  Figures  4.10  and  4.1 1  show  the  results  of  the  flaw  in  the  measurement. 


Figure  4.9:  Poorly  Made  Sample 


To  further  understand  the  sensitivity  to  sample  fabrications,  some  more  specimens 
were  modified  in  exaggerated  ways  to  observe  the  behavior.  The  first  modification  was  to 
cause  a  large  hole  around  the  ridge.  An  otherwise  well  fabricated  sample  was  cut  out  and 
then  modified  as  shown  in  Figure  4.12  and  the  measurements  are  shown  in  Figures  4.13 
and  4.14. 


58 


10  r 


05 

>s 

> 

g 

E 

a) 

Q_ 

<D 

> 

+-> 

TO 

Q) 

a: 


6  - 


4  - 


2  - 


J* 


-  Real  er 
-Imaginary  e 


1 


-I_ I_ I_ I — 


—1_ I_ I_ I_ I_ I 


7  8  9  10  11  12  13  14  15  16  17  18 

Frequency  (GHz) 


Figure  4.10:  Poorly  Fabricated  Sample:  Relative  Permittivity 


Finally,  modifying  the  troughs  (or  ears)  of  the  guide  were  investigated  (Figures  4.15, 
4.16  and  4.17).  Judging  by  the  results,  the  overall  permittivity  was  driven  down  while  the 
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Figure  4.12:  Modified  Ridge  Sample 


Figure  4.13:  Modified  Ridge  Sample:  Relative  Permittivity 


results  in  general  were  completely  unreliable  for  the  higher  frequencies.  This  error  and  the 
associated  uncertainty  could  be  a  limiting  factor  for  making  reliable  measurements  with 
DRWG  particularly  at  higher  portion  of  the  bandwidth. 
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Figure  4.15:  Modified  Trough  Sample 


In  addition  to  the  problems  discussed  above,  there  were  also  challenges  met  with  the 
sample  fitting  too  tightly  in  the  guide.  This  caused  compression  in  the  sample,  which  will 
negatively  affect  the  measurement.  In  addition  to  the  compression  causing  changes  in  the 
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Figure  4.16:  Modified  Trough  Sample:  Relative  Permittivity 


Figure  4.17:  Modified  Trough  Sample:  Relative  Permeability 


constitutive  parameters,  it  was  also  applied  non-uniformly  across  the  sample.  The  results 
were  essentially  useless  since  no  accurate  sample  thickness  could  be  applied. 
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4.5  Comparison  to  Rectangular  Waveguide 

To  validate  the  methodology  presented  here,  the  AFRL  sample  was  also  measured 
with  rectangular  waveguide.  While  the  samples  were  not  the  same  thickness  or  made  at  the 
same  time,  the  results  should  show  similarities.  These  results  are  shown  in  Figures  4.18  and 
4.19.  The  overall  comparison  is  well  within  any  uncertainty.  The  rectangular  waveguide 
seems  to  validate  using  DRWG  as  a  measurement  apparatus. 


Figure  4.18:  DRWG  Compared  with  Rectangular  Waveguide  Measurement  ( er ) 


4.6  Additional  Challenges  in  Measurement 

Unique  challenges  exist  in  DRWG  measurements  with  regard  to  the  higher  measure¬ 
ment  frequencies,  in  comparison  to  rectangular  waveguide.  Initial  results  were  plagued  by  a 
common  decrease  in  the  final  measurement  at  approximately  16  GHz  (Figure  4.20),  which 
is  a  higher  band  than  is  measurable  by  X-band  waveguide.  After  some  extended  experi¬ 
mentation,  it  was  discovered  that  the  uncalibrated  data  displayed  the  same  decrease  (Figure 
4.20),  indicating  a  problem  in  the  measurement  technique  or  a  piece  of  equipment.  Two 
problems  were  identified:  inconsistent  clamping  pressure  and  a  damaged  reflect  standard. 
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Figure  4.19:  DRWG  Compared  with  Rectangular  Waveguide  Measurement  (p.r) 


Upon  close  inspection,  one  of  the  reflect  standards  had  a  slight  curve  in  it  (Figure 
4.21).  This  was  causing  a  slight  cavity  response  and  defeating  the  assumption  of  a  zero 
offset  short  in  the  TRL  calibration.  The  small  curve  was  shown  by  placing  a  razor  blade 
across  the  standard  and  observing  light  coming  though.  The  non-curved  reflect  standard 
shows  a  little  light  coming  through,  where  the  razor  has  some  slight  imperfections.  The 
curved  reflect  standard  shows  a  significant  amount  of  light  passing  under  the  blade,  thus 
illustrating  a  significant  problem. 

While  the  solutions  described  above  seemed  to  mitigate  measurement  problems  in 
the  rigid  samples,  they  did  not  solve  the  all  of  the  problems  with  the  absorbers.  Close 
inspection  of  the  samples  revealed  a  two  major  problems.  The  first  was  simply  poor 
fabrication  that  resulted  in  some  unwanted  gaps  between  the  material  and  sample  holder. 
The  second  problem  was  identified  by  measuring  the  sample  thickness  while  in  the  holder. 
The  thickness  was  actually  slightly  larger  than  when  outside  of  the  guide,  indicating 
compression  of  the  sample.  The  dimensions  of  the  sample  were  reduced  by  carefully 
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Figure  4.20:  Curved  Reflect  Standard  vs.  Flat  Standard  Raw  Data 


(a)  Curved  Reflect  Standard  (b)  Flat  Reflect  Standard 

Figure  4.21:  Reflect  Standard  Comparison 


scraping  material  off  with  a  razor  knife  until  the  thickness  of  the  sample  did  not  change 
when  it  was  inserted  into  the  holder.  Care  was  also  taken  to  take  minimal  material  off  of 
the  ridged  portions,  while  still  slightly  reducing  the  overall  size. 
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V.  Conclusion 


The  method  of  using  DRWG  was  presented  for  gathering  constitutive  parameter  data 
about  a  material.  The  method  required  some  knowledge  about  other  methods, 
particularly  rectangular  waveguide  measurements,  to  appreciate  the  advantages  of  using 
DRWG.  While  a  significant  advantage  of  using  DRWG  is  an  increased  bandwidth 
measurement  range,  the  disadvantages  were  not  trivial.  Material  samples  in  particular 
were  much  more  difficult  to  fabricate.  The  samples  measured  included  a  AFRL  provided 
dielectric  and  an  acrylic  sample  that  both  demonstrated  fairly  constant  constitutive 
parameters  across  the  measurement  band.  The  other  samples  measured  included  two 
commercially  available  microwave  absorbers  that  demonstrated  dynamic  parameters  across 
the  bandwidth. 

In  order  to  understand  the  purpose  of  the  measurements,  a  brief  electromagnetics 
overview  was  presented.  To  appreciate  the  gathering  of  data,  a  calibration  technique  known 
as  TRL  was  discussed  and  derived  using  different  measurement  configurations.  Because 
the  cutoff  frequency  was  important  for  the  NRW  algorithm  to  calculate  the  constitutive 
parameters,  a  detailed  analysis  of  DRWG  was  performed.  Results  were  presented  to  include 
uncertainties  in  the  measurements  and  comparison  to  rectangular  waveguide. 

5.0.1  Sample  Fabrication. 

Sample  preparation  was  one  of  the  most  challenging  aspects  of  the  work  performed 
here.  For  rectangular  waveguide,  samples  can  be  easily  fabricated  from  larger  specimens 
of  the  sample  material  with  straight  razor  blade.  While  AFRL  is  able  to  easily  grow  a 
sample  to  exact  specifications  and  dimensions,  this  may  not  be  practical  in  some  laboratory 
environments  or  possible  with  most  materials.  With  a  soft  material  such  as  the  microwave 
absorber  presented  here,  a  rough  template  can  be  cut  with  a  sharp  blade  and  then  forced 
into  a  line  standard.  A  harder  material,  such  as  the  AFRL  sample,  can  be  formed  in  a 
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similar  fashion,  but  the  template  must  be  more  exact.  Some  tooling  could  possibly  be  used, 
such  as  an  end  mill,  but  the  final  cut  would  have  to  be  forced  into  the  line  standard.  If  this 
method  is  to  be  used,  it  may  be  prudent  to  have  a  line  standard  made  out  of  a  harder  metal 
that  could  prove  to  be  more  robust  when  cutting  samples.  Additionally,  perhaps  some  sort 
of  die  could  be  made  to  cut  the  samples. 

5.0.2  Non  Destructive  Testing. 

One  of  the  major  deficiencies  in  the  type  of  measurement  presented  here  is  the  fact 
that  the  original  test  article  is  destroyed  in  the  process  of  taking  a  measurement.  Harvesting 
samples  from  the  LO  treatment  on  an  aircraft  is  not  usually  done  in  practice,  since  extensive 
and  expensive  repair  work  would  be  required.  In  addition  to  the  need  to  harvest  some 
material,  the  samples  may  not  hold  up  to  repeated  measurements  taking  place.  Because  of 
the  low  tolerances  allowed  around  the  sample  in  the  holder,  a  few  of  the  AFRL  samples 
were  broken  due  to  repeated  extractions  from  the  line  standard.  A  non  destructive  method 
to  measure  samples  could  be  more  convenient  and  would  make  taking  measurements  much 
more  simple  and  has  been  demonstrated  using  rectangular  waveguide. 
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